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The d-dimensional fractional Brownian motion (FBM for short) 
B t = ,B ( t d) ),t £ R) with Hurst exponent a, a € (0, 1), is a 

d- dimensional centered, self-similar Gaussian process with covariance 

E[B^Bl j) ] = !<5i,j(|s| 2Q + I*! 2 " ~ l*~ s | 2 ")- The lon g-standing prob- 
lem of defining a stochastic integration with respect to FBM (and the 
related problem of solving stochastic differential equations driven by 
FBM) has been addressed successfully by several different methods, 
although in each case with a restriction on the range of either dor a. 
The case a = ~ corresponds to the usual stochastic integration with 
respect to Brownian motion, while most computations become singu- 
lar when ct gets under various threshhold values, due to the growing 
irregularity of the trajectories as a — > 0. 

We provide here a new method valid for any d and for a > j by 
constructing an approximation F(e) t , e — » 0, of FBM which allows 
to define iterated integrals, and then applying the geometric rough 
path theory. The approximation relies on the definition of an analytic 
process F z on the cut plane z 6 C \ M of which FBM appears to be 
a boundary value, and allows to understand very precisely the well- 
known (see [5]) but as yet a little mysterious divergence of Levy's 
area for a — » | . 

0. Introduction. The (two-sided) fractional Brownian motion t — ► Bt, 
t G R (FBM for short) with Hurst exponent a, a £ (0,1), defined as the 
centered Gaussian process with covariance 

(0.1) E[B s B t ] = l(\s\ 2a + \t\ 2a -\t-s\ 2a ), 

is a natural generalization in the class of Gaussian processes of the usual 
Brownian motion, in the sense that it exhibits two fundamental proper- 
ties shared with Brownian motion, namely, it has stationary increments, 
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namely E[(B t - B S )(B U - B v )] = E[(B t+a - B s+a )(B u+a - B v+a )} for every 
a, s, t,u,v G R, and it is self-similar, namely 

(0.2) VA>0 (B\ t ,t G R) ^==^ (\ a B t ,t G R). 

One may also define a d-dimensional vector Gaussian process (called: d- 

dimensional fractional Brownian motion) by setting = (B^ , . . . ,B^), 

where [B^ , t G R)i=i,...,<2 are <i independent (scalar) fractional Brownian 
motions. 

Its theoretical interest lies, in particular, in the fact that it is (up to 
normalization) the only Gaussian process satisfying these two properties. 

A standard application of Kolmogorov's theorem shows that FBM has a 
version with (a — e) -Holder paths for every e > 0. In particular, all its paths 
possess finite q- variation for every q > — , in the sense that 



(0.3) sup [2^\B tl -B tl _ 1 \ <l )<oo a.s., 

S = t0< — <t n =t — Q J 

where the sum ranges over all partitions (s = to < t± < ■ ■ ■ < t n = t) of any 
order n of the interval [s,t]. 

There has been a widespread interest during the past ten years in con- 
structing a stochastic integration theory with respect to FBM and solving 
stochastic differential equations driven by FBM. One of the strategies con- 
sists in trying to imitate the Brownian (or, more generally, martingale or 
even semi-martingale) construction by defining, for instance, the so-called 
symmetric integral 

* ,„ , r t x s + X s+£ B s+£ - B s , 

X o dB o = w — lim / ds. 

o e^oJo 2 2e 

This idea and generalizations thereof (by using the classical Newton-Cotes 
approximation scheme for integration) led to the so-called Russo-Vallois 
theory of integration (see, e.g., [9, 18, 19]) which is well defined for the 
scalar fractional Brownian motion with any Hurst exponent a > 0. 

Another approach by using the Malliavin calculus and Skorokhod integra- 
tion is due to Decreusefond and Ustiinel (see [6] ) and Cheridito and Nualart 
(see [4]). Again in the scalar case, they construct an extension of the diver- 
gence operator for every Hurst exponent a > and retrieve, in particular, 
that a = | is a barrier for the existence of the symmetric integral. 

A very different approach, more suitable for nonscalar FBM, is ongoing 
since the work of Lyons; see [14]. This is the one we are going to follow. We 
shall need to recall the main results and notation of [5, 14] before we can 
state our own results. 

Let V be a Euclidean space, and V® k = V <g> ■ ■ ■ (g) V be the kth. tensor 
product of V endowed with a compatible norm || • ||. Lyons starts by defining 
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a space of nonsmooth paths T : M + — ► V on V, called geometric rough paths, 
by completing the space Q°°(V) of smooth paths. For T 6 £l°°(y), k > 1, 
and s, i > 0, we let 



(o.4) r* t = / dt x ® f 1 v t2 dt 2 • • • ® / fc " r; fc dt fc . 

This is an element of with components (r* t)ii,... i* defined in an 

obvious way by projecting r tl , . . . , T' tk onto the corresponding components. 
Repeating the same construction for another smooth path A, one may define 
the following q-variation distance for every q > 1: 

(0.5) 5 q (T,A) Sjt = sup ||A u -r u ||+ sup d k (T,A) Sjt , 

ue[s,t] i<fc<[g] 

where [q] is the entire part of q and d k is the distance defined by 

/ n \ k /l 

(o.6) 4(r,A) M = sup EK_ lit ,-r*_ 1)t ,||* /fc • 

\s=t < --<t n =t l=l ) 

Note that dk is the distance associated to a semi-norm that we shall call 
the dk-norm (which depends on the interval [s, t]). The completion of £l°°(V) 
with respect to the q- variation distance will be denoted by Q q (V). 

The following extension principle is proved in [14]; see also the book by 
Lyons and Qian [15] or the nice review [12]. 

Extension principle. For a given smooth function / :R — > R n — > 
Lin(JR rf ,R n ), let Y be the solution of the stochastic differential equation 
dYt = f(t, Yt) dXf. Then the ltd map 

(0.7) i-.n 00 ^) -^n°°(R n ), x^y 

is continuous with respect to the ^-variation distance. Hence, X admits a 
unique continuous extension X: f2 9 (R d ) — ► f] 9 (IR n ). 

This result (obtained by the usual Picard iteration method) contains, in 
particular, an extension theorem for iterated integrals X — > Jq dX Sl fi (X Sl ) x 
Jo 1 '"Jo" 1 dX Sn f n (X Sn ) . 

Hence, the problem of solving differential equations driven by FBM or 
giving a sense to iterated integrals with respect to FBM may be reduced to 
that of constructing an approximation of FBM which converges in the sense 
of rough paths. 

The paper is organized as follows. We give in the first section a series 
decomposition of the fractional Brownian motion (see Lemma 1.1, Definition 
1.2 and Corollary 1.7), 



(0.8) r ' = E(/ A( s )<k)# + E(j[ M s ) ds )tk~i 



fc>0 v u ' fc>0 
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where £jfc\£fc = £fc are independent complex normal Gaussian random vari- 
ables, and show that Tf may be seen as the restriction to the real axis of the 
real part of a Gaussian process T+ defined on the closed upper half-plane 
n + which is a.s. analytic in z on 1~I + . Equivalently (think of Schwarz's reflec- 
tion lemma) , one may extend r J to the lower half-plane by simply taking its 
conjugate, TJ := (r^~). In this way, FBM appears as the boundary value of a 
Gaussian process defined onC\l. Lemma 1.5 is important for it shows that 
the (a — e) -Holder property of FBM extends to the process on the closed 
upper-half plane. 

Interestingly enough, the T-process may be seen (after a Cayley trans- 
form) as a random entire series living in the unit disk. Lemma 1.5 implies 
then that convergence to FBM on the unit circle is uniform in all directions 
(even tangentially). 

The reader may be puzzled by the fact that we actually obtain a decom- 
position of the fractional noise, that is, the (distribution-valued) derivative 
of FBM (hence the integrals over s in the above definition of Tj). This is 
directly related to the stationarity (and simplicity of course!) of the "covari- 
ance kernel" which is (up to a coefficient) the function (x,y) — > \x — y\ 2a ~ 2 
(at least for a > \ , for otherwise it is singular) . This kernel has many known 
decompositions in terms of orthogonal polynomials or special functions; the 
choice we made for the functions fk is the simplest one, but there are other 
natural possibilities which also lead to a straightforward analytic extension 
to the upper half-plane. This approach may be explained in "physical" terms 
by saying that we approximate the singular kernel (x — y) 2a ~ 2 by a regu- 
lar kernel ((x — y) ± ie) 2a ~ 2 — or, more precisely, (±i(x — y) + e) 2a ~ 2 (see 
below) — where e is a short-distance cut-off. 

The paper is organized as follows. 

Section 1 is devoted to the construction of the T-process, defined as the 
above series (0.8), and to the proof that its boundary value on M. has the 
same law as FBM. The main result of the section is Lemma 1.5. 

Section 2 can be skipped for a first reading since the results will not 
be used in the following sections. It is devoted to the proof of convergence 
theorems for the series. It may be of interest for those interested in simulating 
FBM. Unfortunately, it is absolutely of no use when trying to construct a 
stochastic integration with respect to FBM, except maybe in the easy case 
a > 2 where simpler methods are more adapted. 

In Section 3 we introduce our e-approximation T(e)t of FBM. The process 
is simply defined as the real part of Tf, ie . The e-shift of the variable t into 
the upper half-plane transforms FBM into an analytic process. Contrary 
to the usual schemes of linear interpolation, exact computations are still 
possible when replacing FBM with r(e)j. 

Section 4 is the hard part. It is devoted to the proof of the existence 
of the rough path limit of T(e) for e — > when a > j, which implies the 
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possibility to define stochastic integration with respect to FBM and to solve 
FBM-driven stochastic differential equations. It is mostly interesting from 
a practical point of view, since the existence of a rough path limit for a > 
| is already known from the paper by Coutin and Qian [5]. Namely, we 
compute explicitly the second moment of the Levy area of r(e) in terms of 
hypergeometric functions (see Theorem 4.4), from which the divergent terms 
in the limit e — ► (for a < j) may be precisely identified. We also compute 
the second moment of the Levy volume of T(e) for a three-dimensional FBM. 
It is less explicit, given as a sum of terms, each term being an integrated 
hypergeometric function with respect to a power kernel (see Section 4.3) 
which converges either for a > j or for a > g. 

It would be interesting to know whether the piecewise linear approxima- 
tion used in [5] and the analytic approximation defined in this paper define 
in the limit the same stochastic integration theory and give the same solu- 
tion to stochastic differential equations. We shall come back to this at the 
end of Section 4.1. 

Note finally that Feyel and de la Pradelle (see [8]) proved an ltd formula 
for FBM (actually, for the Liouville fractional Brownian motion, defined as 
fractional integrated Brownian motion) by using an analytic continuation in 
the parameter a for Re a > \. Their approach is totally different, so are the 
results. 

Notation. We will use the following notation throughout the article. 

The Hurst exponent a 6 (0, 1) is fixed once and for all. It is intended 
implicitly to differ from \,\,\-,--- unless explicitly mentioned. One could 
actually suppose that a < \ since there are no singularities at all for a > \ 
(and, hence, the procedure of analytic extension to the upper half-plane 
is pointless), and the theory of stochastic integration of Brownian motion 
(corresponding to the case a = ^) is pretty well understood without all this 
machinery. Nevertheless, all computations below also hold true for a > | 
(unless explicitly stated). 

The one-dimensional fractional Brownian motion process with Hurst ex- 
ponent a will be denoted by (Bt,t £R). A d-dimensional FBM with inde- 
pendent exponents will be written as B = (B^\ . . . ,B^). Alternatively, a 
process with the same law as (Bt , t £ R) may be written as (r^ , t G R) if nec- 
essary. The process T will be obtained as the boundary value of an analytic 
process defined on C \ R. 

If x 6 R and k = 0, 1 . . . , then (x) k = x(x + 1) • • • (x + k - 1) = is 
the Pochhammer symbol. 

If x £ M, then E(x) = sup{n E Z|n < x} is the entire part of x. 

If z G C\ R_, and f3 G C, then z@ := exp/?ln(z) is defined by using the 
usual determination of the logarithm. 
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If x G R, x ^ 0, then sgn(x) = A is the sign of x. 

The Poincare upper half-plane is denoted by II + = {x + iy G C|?/ > 0}, its 
closure by fl + = {x + G <C|y > 0}. Similarly, II - and ft - denote the open, 
respectively, closed lower half-planes. 

In the next proposition, we summarize for the unfamiliar reader the few 
things that should be known on the Cayley transform. 



Proposition 0.1. The Cayley transform C :t — > z = |r| is a bijection: 

(i) from the one-point compactification K of the real line onto the unit 
circle; 

(ii) from the set {t G C| Imt > 0}II{oo} (equal to the one-point compact- 
ification ofti + ) onto the closed unit disk D. 

In particular, C(0) = —1 and C(oo) = 1. The inverse transform reads t = 



1. Definition of the T-process and connection with the first chaos of FBM. 

For a > i , the first chaos of FBM may be constructed by defining TL a to be 
the Hilbert space completion of the space of step functions with respect to 
the scalar product 



(1.1) (f,g) Ha = a (2a-l) I I f(x)g(y)\x-y\ 2a - 2 dxdy 
and setting 



E 



f{t)dB^j (f g{s)dB i 



if, 9) 



Ho 



for /, g G TC a . The space 7i a may be called a reproducing kernel Hilbert space 
for FBM. Note that the convolution kernel K(x,y) = \x — y\ 2a ~ 2 is in Lj oc 
only in the range a > |. If a < s> then this definition does not make sense 
any more and the construction of a reproducing kernel Hilbert space is more 
complicated, relying, for instance, on a mapping to L 2 (W) (and hence to 
Brownian motion) by means of fractional derivatives; see [17]. 

One would like to define the "derivative" of FBM as "the" process with 
covariance E[S^,S^] = a(2a — l)\x — y\ 2a ~ 2 , but this has no meaning at all 
because of the nonintegrability of the singularity on the diagonal x = y. Note 
also the awkwardness of the absolute value in the kernel \x — y\ 2a ~ 2 , which 
is, however, necessary since the power functions are multivalued. 

By going over to the complex plane, all these problems may be avoided, 
as we shall show presently. 
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Lemma 1.1. Let fk,k > be the analytic functions defined on the upper 
half-plane as 



h{z) = 2 a-i . Ml -2a) J (2- 2a) k (z + i\ 2a ~ 2 (z - i 



2cos7ra V k\ \ 2i I \z + i 



If z,w G n+ ; then the series J2k>o fk( z ) fk(w) converges in absolute value 
and 

fk{z)7kiw) = ^—^-(-i(z - w)f a ~ 2 . 

~ n l COS TTCt 

Remark. The above definition makes sense since Re(— %{z — w)) > 
(recall from the Introduction that the fractional powers are defined by means 
of the usual determination of the logarithm on C\ R_ ) . 

Proof of Lemma 1.1. By Proposition 0.1 in the Introduction, the 
Cayley transform z — ► |t4 maps the upper half-plane onto the unit disk, 
hence, the convergence of the series S(z,w) := ^"f^a) J2k>ofk( z )fk(w). An 
explicit computation yields 

z + i\ 2a - 2 f w - i\ 2a ~ 2 f z-iw + i^ 2a - 2 



1-2) p-I p- 1 



2i J \ —2i J \ z + iw — i 
= (-i(z-w)) 2a ~ 2 . □ 

Definition 1.2. Let A; = 0,1,... be independent centered real 

Gaussian random variables with variance 1, and let £ fe = ^ ± t£|. The series 

(z g n + ) 

fc>0 

defines a random process r'_, + ; we denote by T'f its complex conjugate 
Efc>o /fcWfjfe • 

With this notation, r' + (resp. T'~) lives on II + (resp. II - ). The above 
lemma yields immediately the covariance 

E[r ,+ (z)r ,+ (w)] = E[r'- (z)r'- («;)] = o, 
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By using the Cayley transform z — > C = §rf> the process T /+ may be 
rewritten (up to a prefactor) as a random entire series 1 , namely, 




(2-2a) fc fc + 
2cos7ra fc! ^ ?fc 

Random series have been extensively studied; see, for instance, the book 
[10] by Kahane. The facts we need on random entire series are elementary, 
so we chose to give self-contained proofs. The next lemma will allow us to 
prove that r /+ is analytic in the upper half-plane. 

Lemma 1.3. 

1. Let (a n )n=o,i,... be independent random variables. Then there exists R> 
such that f(z) = ^2 n >oZ n a n (z 6 C) converges a.s. for \z\ < R and 
diverges a.s. for \z\ > R. The function z — > f(z) is a.s. analytic for \z\ < R 
and the convergence is uniform on every compact C B(0,R). 

2. The radius of convergence R is given by 



^ = sup(A>o|^P[|a n |>A re ] = ool. 

I n>0 J 



Proof. 



1. By Hadamard's criterion, the radius of convergence R of a series J2n>o a nZ r ' 
is given by 

1 



limsup \l\a T 

ti n 



Kolmogorov's (0-1) law states that P[limsup ^/|a n | > A] = or 1 since 
this is a tail event. Hence, l/R = sup{A > 0|P[limsup \/|a n | > ^4] = 1} if 
this set is nonempty, or else. 
2. By the Borel-Cantelli lemma, 



PUimsup sj\a n \ > A] 



n>0k>n 



1 



if and only if En>o p [KI > An ] = °°- □ 

Corollary 1.4. Suppose f(z) = J2 n >o ^nZ n £, n , where the (£ n ) are inde- 
pendent standard Gaussian variables. Let l/R := limsup n ^ oc \/|A n |. Then 
the series converges uniformly on every compact C 13(0, R) to a holomorphic 
function. 



iOne should actually set f'+(C) = F'+(z)^ = ^rfr'tO) since r ' + is to be in- 
terpreted as a derivative process. But this is not important for what follows. 
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PROOF. Supposing A > l/R, then P[|A n £ n | > A n ] < C( e A n/ Xn " ) goes 
rapidly to and the series converges, while it is clear that X)n^[l^n£™| — 
A n ] =00 if A< l/R. □ 

It follows immediately from the preceding corollary and (1.3) that T'+ is 
well defined and analytic for z S II + . 

The next technical lemma is crucial to establish the connection with FBM 
and regularity properties. 

Lemma 1.5. Let 7, 7': (0,1) — > n + be continuous paths with endpoints 
7(0) = a\ + £61,7(1) = 02 + ib 2 , 7'(0) = ct'i + i6' 1 ,7 / (l) = a' 2 + ib' 2 in the closed 
upper half-plane tl + , that is, 61,62,6^,62 > 0. Then, for every a E (0,1), 



1. T/ie double integral 



I = I du l_dv(—i(u — v)) 2a 2 



is well defined. It depends only on the homology class of the paths 7,7', 
or, in other words, on the endpoints 7(0),7'(0),7(1),7'(1). Furthermore, 
it is invariant under real translations 7 ^ 7 + a, 7' — > 7' + a (a€M.). 

2. Suppose that a\ = a[ = and take 61, 62, b[, b' 2 — ► 0. T/ien 

e i7rasgn(a' 2 -a 2 )| a ^ _ a2 | 2Q _ g- i7ras g n ( a 2) | a2 |2a _ giTrasgr^aj) | 1 2a 

/7 ' 7 ' ^ 2a(2a - 1) ' 

Hence, 

- — 2cos7ra |a 2 | 2a + |a 2 | 2a - \a' 2 -a 2 \ 2a 

7 7j y +j 7j y-2Ke7 7 ,y ^ a(i _ 2a) • 2 

3. Suppose 7 = 7'. T/ieri i/iere exists a constant C > suc/i i/iaf 
du / - v)) 2a ~ 2 < C| 7 (l) - 7 (0)| 2a . 

Proof. Point 1 follows from Cauchy's theorem and the fact that 
\{-i{u - v)) 2a ~ 2 \ <{lmu + lmv) 2a - 2 

and 

ete / dy(x + y) 2a ~ 2 < 00. 
Jo 

As for point 2, an explicit computation yields (if a\ = a± = 0) 
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_(— i(u — v)) 2a 2 dudv 



'7 J f 

= ((6a + b' 2 - i{a 2 - a' 2 )f a - (6 2 + b[ - ia 2 ) 2a 

- (6i + b' 2 + m' 2 ) 2a + (&! + 6' 1 ) 2Q )(2a(2a - l))" 1 . 
When b±,b 2 , b' x ,b' 2 — ► 0, this expression goes to 

e i™sgn(4-a 2 )|^ _ a2 |2a _ g-i7rasgn(a2) | fl2 |2a _ e j7rasgn(a 2 ) | a / 1 2a 

2a(2a - 1) ' 

hence, also the result for Rel 7i y. 

So let us now set about to show point 3. We split the proof into several 
parts. 

1. Suppose 7(0) = and 7(1) = z, Imz > 0. Then one finds 

dudv(-i(u - v)) 2a ~ 2 = (2lmz) 2a - 2Re(-iz) 2a 

= {2rsm6) 2a - 2Rer 2a e 2ia ^- n/2) , 
if z = re ld , 9 S [0, it], hence, there exists a constant C > such that 



(1.4) 



7 J i 



\2a-2 



< C\z 



2a 



dudv(—i(u — v)Y 

By invariance under real translations z — > z + a, a E R, a similar estimate 
holds for 7(0) real. 

2. Suppose ^2 < ^1 < ^2- Set it = xi + iy\, v = Xi + iy2- Let us compute the 
above double integral I 7i7 by choosing 7 = 71 U 72 := [a\ + ib\,a 2 + i&i] U 
[02 + ib\,a 2 + ^2] as a horizontal line followed by a vertical line. 

The double integral on the vertical line / L 2X72 dudv(—i(u — v)) 2a ~ 2 
is less than 

fe 2 /-62 

/ dy 1 dy 2 (2b 1 ) 2a - 2 

2q-2/, , n2 /l^2 - &lh 2_2a ,, , 1 2a 



(1.5) <C6f-^6 2 - ^=1^-^) I62 — 61 



< |&2-&i| 2Q . 



The mixed integrals 



Jy := / / dudv(-i(u-v)) 2a 2 , ij^j, 
7i •'Tj 
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satisfy 

r\a.2—ai\ rb2~b\ 

\hj\<C dx dy(2bi + x + y) 2a ~ 2 

Jo Jo 

^C^-a^-b^b 2 *" 2 



(1.6) 

hence, also 
(1.7) \Iij\<C'\a x -a 2 

On the other hand, 

Q2— ail rb2—bi 



. 2a-l 

ai — a 2 



r\a,2— ai /•02—01 

/ / dy(26i +x + y) 

jo Jo 



2a-2 



(1-8) = _ 1} [(61 + 6 2 + |a 2 - ai |) 2a - (&x + 6 2 ) 2 " 

-(2b 1 + |a 2 -a 1 |) 2 - + (26 1 ) 2 "]. 

We must now distinguish 4 cases according to the sign of a — 5 and to 
the relative order of magnitude of b\ and |ai — a 2 1 . 

Suppose first a > ^: then < C'\a\ — a 2 \ 2a by (1.7) if 61 < \a\ — a 2 | , 
and \Iij\ < C'\ai - a 2 \ 2a ~ 1 (b 2 - h) if h > \ ai - a 2 \ by (1.6). Holder's 
inequality 

\x\p \y\i 1 1 

\xy\ < 1 , - + - = 1 

p q p q 

applied to x = \m - a 2 \ {2a ~ 1)/a , y = {b 2 - 6i) 1/a , p = ^1 yields in the 
latter case |jy < C"(\a x - a 2 \ 2 + |6 a - 6 2 | 2 ) a . 

If now a < i, then |^-| < C>i - a 2 | 2a by (1.7) if is bounded 

below by a positive constant, while (1.8) yields for |&i| <C \a± — a 2 \ 

< Csup(6 2Q , |ai - a/ - 1 ^ - 61)) < |ai - a 2 | 2c \ 

Finally, the integral on the horizontal line 



hi = I I dudv(—i(u — v)) za 1 = I I dxidx 2 (2bi — i(xi — x 2 )) 2a 2 

J71 •'71 J ai J a\ 

can be evaluated explicitly, namely, 

= 2(2b!) 2a - (2b! - i(oi - a 2 )) 2a - (2&i + ijax - a 2 )) 2a 
11 ~ 2a(2a - 1) 

An expansion at order two by using Taylor's formula gives In < Cb 2a ~ 2 \ai — 
a 2 | 2 < \a± — a 2 \ 2a if b\ \ai — a 2 \. If 61 is of the same order as a± — a 2 or 
smaller, the same estimate comes out trivially. 
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3. Suppose now that b 2 > 2&i. Choose the following contour of integration: 
set 7 = 71 II 72, with 71 = {a\ +ib\{l — t)\t G [0, 1]} and 72 = {(1 - t)ai + 
t(a 2 + ib 2 )\t£ [0,1]}. Then 

dudv(-i(u - v)) 2a ~ 2 = In + I12 + hi + ^22, 

In/ J-y 

where Ijk = / 7 . I lk dudv(—i(u — v)) 2a ~ 2 . The integral In (resp. ^22) is 
less that a constant times [resp. {(a\ — a 2 ) 2 + by point 1. As for 
the mixed terms Iij,i 7^ j, an explicit computation yields, for instance, 

h2 = rr[(6i + 6 2 - i(ai - a 2 )) 2a - - (6 2 - i(ai - a2)) 2 °], 

2a{2a — 1) 

which is < C(\ai — a 2 \ 2 + \b\ — b 2 \ 2 ) a . By putting everything together, 



one gets 



dudv(—i(u — v)) 2a 2 



<C'((ai-a 2 ) 2 + (bi-b 2 ) 2 ) 



2\a 



□ 



Definition 1.6. Define for any z e II+ 

duT'+, 



where 7 : (0, 1) — ► II + is any continuous path with endpoints 7(0) = 0, 7(1) = 
z and, similarly, TJ = (Tf). 

Lemma 1.5 shows that has a well defined covariance on YL^ 1 . One has 
thus obtained the r ± -processes. These are hermitian Gaussian processes 
defined on TL^. Restricted to H^, they have analytic paths. 

Corollary 1.7. 

1. The real-time real-valued process (Tt,t > 0) defined as 

r t = 2Rer+ = r+ + r7 

has the same law as (B t ,t £ R). 

2. Let t € R. Then, with probability 1, 

2R e r+ = r+ + rj ->r t 

as z € II + goes to t along any path in the closed upper half-plane. More 
precisely, 

E[\ReTt -T t \ 2 ]<C\z-t\ 2a . 

Hence, the real-valued process z — ► , z £ II + has an (a — e) -Holder 
continuous version for every e > 0. 
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Remark. In other words, FBM appears as the boundary value of the 
analytic process defined as on Il + and Tj on II - . We shall come back 
to this in Section 3. 

Proof. It is straightforward from Lemma 1.5 by using Kolmogorov's 
lemma for Gaussian processes. □ 

We shall see in Section 3 that the convergence of Re Tf to FBM as Im z — > 
is a.s. uniform on every compact. 

But let us first come back to the series T'+ = J2k>o fk(z)(, k i z £ n + ) and 
show how FBM appears as the limit of a series defined directly on M. 

2. On a series decomposition of FBM. Let us recall the definition of 
fk(z),z G II + (k > 0) from Lemma 1.1 in Section 1: 



fk(z) 



■><*— l . 



2cos7ra V k\ \ 2i I \z + i 



'a(l-2a) 1(2 -2a) k fz + i\ 2a ~ 2 f z - i x k 



This function extends analytically to the half-plane Imz > — 1; in any 
case, it is well defined on the real axis. By Stirling's formula, 



(2-2a) k _ I 1 r(2-2a + fc) _ 1/2 

V k\ -\/r(2-2a)V r(i + fc) -° k ^ k >■ 

Recall also that = 1 for z G R. Hence, the series J2k>o |/fc(*)| 2 diverges 
on the real axis as X}fc>i^ 1_2Q '■> which is of course not at all surprising 

[otherwise the random series J2k>o fk(t)it + Sfc>o /fc(*)£fc would have been 
a good candidate for a would-be derivative of FBM!]. 
Set 



F k (z):= f f k (s)ds, z G II + . 
Jo 



Lemma 2.1. 



1. Let C a (s,t) = J2k>o Fk( s )Fk(t), s,t 6l. Then this series converges in 
absolute value. 

-zirasgn(s)| „|2a , lira sgn(t) U| 2a _ lira sgn(t-s) u_ _|3a 

2. One has C a (s,t) = - [£J —^ 2 ^ - IL^L_. 

u V ' / 4 cos 7ra 

3. Let B t :=Efc>o^ 1 fc(*)Cfc" + I]fc>o^fc( t )?fe • r/ien (5 t ,tGM) is a two-sided 
FBM with Hurst index a. 

Proof. Let z stay in a fixed compact set K C II + . Then 

,'u — i\ k f u + i s 

du 



o \u + i J \ 2i 
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1 f z , d 
du- 



k + 1 Jo du 



i\ k+1 ~\ (u + i) 2a 
u + i 



(20 



,"\2a-l 



-(2iV- 2a \ — ( z + i) 2a - (_i)fe+i e ^ 
fe + 1 I \z + V 



Hence, 



which guarantees that the series C a (z,w) := J2k>o Fk{z)Fk(w) is normally 
convergent in K. Consequently, 

C a (s,t)= lim C a (z,w), 

z— >s,w— *t 

if z — ► s,w — ► i while staying inside n + . But (see Lemma 1.2) the series 
J2 fk(z)fk(w) converges in absolute value if z, w € n + , hence (by commuting 
the sum with the integral), 

q(l-2a) r* 2Q _ 2 

C Q (2:,u;J = du dv{—i(u — v)) 

2 cos 7TQ; Jo Jo 

Now one may conclude by using Lemma 1.5(2). □ 

By standard arguments using martingale inequlities and Fernique's lemma 
(see, e.g., Theorem 3.3.2 in [2]), the Karhunen— Loeve type series 

(2.i) *^£^)4 + + E*U*)£fc 

fc>0 k>0 

converges uniformly with probability one to t — > Bt on every compact set 
KcM.. 

Let us now study the rate of convergence of the above series. It may be 
computed by using standard entropy arguments that we reproduce from [2] . 
Set i?jv(i) := J2k>N -^fc + J2k>N Fk(t)£k- Then the standard deviation 
semi-metric 



d N (s,t) := ^E[(R N (t) - R N (s)Y 



defines a natural metric on R from the point of view of the Gaussian process 
t — > Risr(t). We shall be using the following 
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Proposition 2.2 (See [20], page 101). Let X t be a centered Gaus- 
sian process on K C M compact, and d(s,t) := ^/E[| — X s \ 2 ] its associated 
(semi-) metric. Then 



E 


sup \X t \ 






.t&K 


Jo 



»+oo 

\nN{d,K,e)de, 

where C is a universal constant, and N(d,K,e) is the smallest number of 
balls of d-radius e covering K. 

The logarithm of N(d, K, e) is usually called the entropy of the process 
X on K. 

Let s,t remain in a fixed compact K C M. Then \R^(t) — Rn(s)\ may be 
estimated either by remarking that 

E[\R N (t) - R N (s)\ 2 ] = 2 ]T \F k (t) - F k (s)\ 2 < 2 £ \F k (t) - F k (s)\ 2 

k>N k>0 

(2.2) 

= E[\B t -B s \ 2 ] = \t-s\ 2a 



or by using an upper bound appearing in the proof of Lemma 2.1, namely, 
E[\R N (t) - R N (s)\ 2 } < 8 £ sup(|F fe (i)| 2 , \F k (s)\ 2 ) 

k>N 

(2.3) 

<C k' l ~ 2a = C'N~ 2a . 

k>N 

Hence (the following trick is borrowed from a similar proof in [7] ) , for every 
couple (p, q) of positive numbers, 

(2.4) E[\R N (t) - R N (s)\ 2 ] p+q < C"N- 2pa \t - s\ 2qa . 

Consequently, lnN(d,K,e) ~ l n (e~^ p+q ^ qa N~ p ^ q ) for e < N~ p a /(P+i) , an d 
In N(d B ,K,e) = for e > JV-WCp+s), leading to 



E 



SUp i?Ar(t) 

teK ' . 



<CjP+AN- pa / p+q [ t/ln(l/x) dx<C'J l±lN- pa / p+q . 
y qa Jo v V qa 



The best estimate for sup tg ^- \R]\[(t)\ is obtained in the limit p — > oo (with 
q fixed), namely 



E 



SUp|i?Ar(i)| 

teK 



for every e > 0, with a constant C e that diverges when e — > 0. 
Summarizing, we have proved the following: 
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Theorem 2.3. Let K C M 6e a compact. Then, for every e > 0, 



TV 



E 



sup 



fc<7V 



fc<Af 



iV— »oo , 



Another unrelated series expansion of the fractional Brownian motion 
(using zeros of Bessel functions) has been obtained by Dzhaparidze and Van 
Zanten in 2004 (see [7]). The convergence appeared to be optimal in the 
sense of a paper by Kiihn and Linde (see [11]): they proved that no series 
converging uniformly on every compact to FBM a could have a remainder 
of order less than N~ a ^\og N. Hence, the convergence is considered to be 
optimal if the remainder is < C £ N^~ a ~ £ ^ for every e > 0. We refer to these 
two papers for details. Here also, as we have just shown, the convergence is 
optimal. 

3. First approximation of FBM. From the results of Section 1, it is nat- 
ural to approximate FBM by the following process: 

Definition 3.1. For any e > 0, let (T(e) t ,t e R) be the Gaussian process 
defined by 

r(e) t = ReT + (t + ie) = T + (t + ie) + T~(t - ie). 

By Lemma 1.5, T(e)t has a (Gaussian) limit Tf in 1? as e — > 0, and the 
limit process t — > Tt has the same law as FBM. In other words, viewing r + 
and r~ as the expression on n + and on n~ of the same analytic process 
r(z) defined onC\l, FBM may be seen as the boundary value of T, namely, 



(3.1) 



T t := lim(r(t + ie) +T(t-ie)). 

e^0 



More precisely, still by Lemma 1.5, there exists a constant C > such that, 
for every t 6 M, 

E[\r(e) t -r t \ 2 ]<Ce 2a . 

3.1. Uniform approximation of FBM by T(e). We shall give in this para- 
graph, by a method similar to that of the preceding section, an estimate of 
the sup- norm of the difference T — T(e). 



Theorem 3.2. Let K C M be a compact. Then, for every n > 0, 



jq-a 



E 



sup|r t - r(e) t \ 

K 



0. 



A STOCHASTIC CALCULUS BY ANALYTIC EXTENSION 17 

Proof. The proof goes as in Theorem 2.3. Namely, the estimate (2.2) 
may be replaced by 

E[\(T t -T(e) t )-(T s -T(e) s )\ 2 ] 

< 4su P (E[|r t - r s | 2 ],E[|r( e ) t - r(e) s \ 2 }) 

<C\t-s\ 2a 
(see Lemma 1.5) and the estimate (2.3) by 

E[\(T t -T(e) t )-(r s -r(e) s )\ 2 ] 

< 4sup(E[|r t - r( e ) 4 | 2 ],E[|r s - r(e) s | 2 ]) 

< C'e 2a 
(also by Lemma 1.5). □ 

3.2. Convergence of analytic iterated integrals. Obtaining good estimates 
of iterated integrals is the key point in order to apply rough path theory. 
We claim that the analytic iterated integrals, as defined below, converge in 
L 2 (P) for every Hurst index a > 0. 

Definition 3.3. Let s, t > and fx, . . . , f n (n > 1) be n analytic func- 
tions defined on a neighborhood of the closed strip = {z = x + iy\0 < x < 
t,0 <y < s}. Then an analytic iterated integral is an integral of the form 

/ d Xl f( Xl + iei)r , (ei)+'(") f * dx 2 
Jo Jo 



(3.2) 



dx n f{x n + i£ n )r ( y £ n ') x ^ n \ 



where s±, . . . ,e n > and 1 < i\, . . . ,i n < d (recall d is the number of compo- 
nents of FBM). 

Theorem 3.4. Let s, t > and fi, ■ ■ ■ ,f n (resp. g\, . . . ,g n ) be n analytic 
(resp. anti-analytic) functions defined on a neighborhood of the closed strip 
f2 = {z = x + iy\Q < x < t, < y < s} (resp. £1). Then 

V(e,rj) :=E Qf dx^xi + ie^fa)^ j*' 

(3.3) 



dx 



2 



x ( / dy igi ( yi -iei)r'(? ?1 ) / dy 



.VI 

S/n-l 



2 



dx n g n (y n - ir] n )r ' {r] n ) y ^ n ' ) 
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admits a limit V(0, 0) when e = (ei, . . . , £ n ),V = (wi, ■ ■ ■ , Vn) go to zero, which 
is bounded by 

n n 

|V(0, 0)1 < CY[sup\fi(z)\ J] sup \ 9i (w)\ 

(3.4) 

x max(s 2Q , ts 2a ~\ t 2a ~ 2 s 2 , t 2 s 2a ~ 2 ) n 
for a universal constant C. 

Proof. Supposing i\ ^ i 2 7^ • • • 7^ i n , one may set j\ = a(ii), . . . ,j n = 
a(i n ) for a certain permutation a [otherwise V(e,r/) = 0]. Let s' = s — sup k e k — 
sup fc r/ fc . Then 

fail — 2a) \ n ( /"* f Xl rxn-i 

V(e,rj) = [-p. / dxi dx 2 --- dx n 

V 2cos7ra / \Jn Jn Jo 



(3.5) 



ryi rvn-i 
x ( / dyi I dy 2 --- dy n 



\Y[f(x k + ie k )j [Y[g(yk-im)j 

~\{-i{x k - y a{k) ) +e k + rj a{k) ) 2a ~ 2 



x 

k 



may be extended as a holomorphic function of x\, . . . ,x n on the strip {z = 
x + iy\0 < x < t,0 < y < s — sup k e k }, respectively, as an anti-holomorphic 
function of y±, . . . , y n on the strip {z = x — iy\0 < x < t,0 < y < s — sup k r] k }- 
Set 

(3.6) jo = [0,is'}U[is',t + is'}U[t + is',t], 7e = 7o + is 

and, if z £ j E , let j e (z) be the section of the path 7e comprised between 
and z. Then, by a deformation of contour, one gets 



VM) 



a(l-2a) 

2 COS IT a 



(3.7) 



dz\ \ dz 2 --- dz n 

lex J7 ei (zi)+i(e 2 -£l) ^7e n _i(^n-l)+*(£n-en-l) 



dw\ I dw2 • • • / dw n 

TVT "'7»Jl(«'l)+i(rp— m) •'7?7 7 i_l( lu n-l)+«(»?ri-»?n-l) 



~[[f(zk) ~[[g(wk) n(-i(zk - w a{k) )) 



2a~2 
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Let M :=Uk sn Pzen\fk(z)\Uk su Pwen\9k(w)\. Then 

|V(e,77)|<C-M(7 \d Zl \--- [ \dz n \]([\dw 1 \---[_\dw ri 



(3.8) 



'7o ^70 / \ J lo J lo 

x ni(— *(^*: -w a{k) )f a ~ 2 \ 



k 

The different components decouple completely and one is left with the prob- 
lem of estimating 

(3.9) I:= [ [ \dz\\dw\\(-i{z - w)) 2a ~ 2 \. 

J 70 J 70 

Splitting 70 and 70 into their three linear pieces, one sees that 



(3.10) 



\dz\ I \dw\\{-i{z-w)) 2a - 2 \ = / dx / dy{x + y) 2a - 2 

[0,is] J[0,-is] JO JO 



n2a o 

„2q. 



2a(2a - 1) 



,2a-2< ^ . 2a-l. 



(3.11) / \dz\ \dw\\{-i{z - w)Y a - z \ <ts< 

J[is,is+t] J[0,—is] 

\dz\ f \dw\\(-i(z-w)) 2a - 2 \<t 2a - 2 s 2 ; 

(3.12) J[o,is] J[-is+t,t] 



(3.13) 



/ \dz\ 

J [ts,is+t\ 



\dw\\(-i(z-w)) 2a ~ 2 \<2 2a - 2 t 2 s 2a - 2 , 

[—is,—is+t] 



hence the result. 

If now some of the indices i\, . . . ,i n coincide, then V(e, rf) may be obtained 
as a finite linear combination of terms of the same type as (3.5), hence, the 
same estimates hold. (Easy) details are left to the reader. □ 

In particular, the second moment of any analytic iterated integral con- 
verges when e — > 0. Hence the divergence of Levy's area for FBM with Hurst 
index a > j stems from the mixed integral Jq dxiT'(ei) + '^ J Xl dx2^' '(£2)" "'^* 2 ' 'j 
see computation of the integral V2 in the proof of Theorem 4.4 below. 

4. Existence of the rough path limit for a > ^. The purpose of this 
section is to prove that T(£) ((i-dimensional £-regularized FBM) has a rough 
path limits when e — ► for every d > 1 and a > \. Let us be a little more 
specific about what has to be done: 
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- for every a > 0, di(T(e), T) s j — > with probability one (which is a simple 
consequence of Lemma 1.5); 

- if a > i, then no regularization is needed. Hence, we assume (for the 
proof of the rough path convergence) that a£(^). Yet the explicit com- 
putation (see Section 4.1) of the e-regularized Levy area for FBM is valid 
for every Hurst index a; 

- if | < a < one needs to prove that T(e) converges in the ^-distance 
for q > 1/a. Since a simple integration by parts yields Y{e) 2 i = \(T(e)^) 2 , 
the only real problem is to prove the convergence of T(e) 2 j for i ^ j. This 
we do in all details in Section 4.2; we shall actually prove that convergence 
in the ^-distance holds for every a > \\ 

- if ^ < a < 3, there remains then to prove the convergence of the third 
iterated integrals r(e)? • k in the e^-distance. Since T(e)f Li = |(r(e)W) 3 and 
r(e)? ij -, T(e)f j i , r(e)| )i)i (i / j) reduce easily (by integration by parts) to 
second-order iterated integrals of one of the following types, 

Jo 

JO 

Jo 

we shall only need to consider the case j ^ k. 

4.1. Computation of two-dimensional Levy's area for a £ (0, 1). Gener- 
ally speaking, given a two-component stochastic process Xt = , -X* ) 
with regular enough trajectories (say, for smooth trajectories), we call Levy's 
area 2 of X between time s and time t the random quantity 

Axe as , t (X)=j\x®-XP)dxU 
(4.1) =j\x^-X^)±X^du 

= f X'Pdx \ X xfUy. 

J s J s 

2 The usual definition of Levy's area is antisymmetric in X^ and X*- 2 ' , given by 
/* Xi 2) dXi 1] - dX™ . It differs from our definition in an inessential way for our 
purposes. 
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It is given geometrically as the area comprised between the curve 




and the union of two broken lines, respectively, horizontal and vertical, 



The purpose of this paragraph is to give an explicit expression of this 
quantity for X = T(e), where T is a two-component FBM. 

We shall be using a number of times the following integral representation 
of Gauss' hypergeometric function 2F1 (see [1], (15.3.1)): 



valid if Rec > Re b > and z € C \ [1; +00). Recall that 2-^1 (« ; b, c; z) (c 7^ 
0, — 1, . . .) is defined around z = by an infinite series with radius of con- 
vergence 1 and has an analytic extension to the cut plane C \ [1; +00 [. The 
connection formulas give 2-^1 ( a > b, c; z) in terms of a linear combination of 
hypergeometric functions in the transformed argument <j)(z), where 4> is any 
projective transformation of the Riemann sphere preserving the set of singu- 
larities of the hypergeometric differential equation, namely, {0,1, 00}. They 
relate the behavior of the hypergeometric functions around with their be- 
havior around 1 and 00 . We reproduce here for the convenience of the reader 
two connection formulas, relating the behavior around with the behavior 
around 00, respectively, 1 (see [1], Section 15.3): 




(4.2) 




2 Fi(a,b,c;z) 



(4.3) 





2 F 1 (a,b,c; z) 



(4.4) 



T{c)T{c-a-b) 
T(c-a)T(c-b) 



2i*i(a, 6; a + b — c+ 1; 1 — z) 




z [l;+oo). 
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Let us also recall that 2-^1(0, b, c; z) is symmetric in the arguments a, b, 
constant (equal to 1) if a = or b = 0, and that lmx,^ ljS ^[ 1;+00 [2-Fi(a, b, c; s) 
exists and is finite if and only if Re(c — a — b) > 0, in which case 

(4.5) 2 Fi(a,&, C ;l) = ^ c)r(c ~ a ~^ , c <£Z_, Re(c - a - b) > 0. 

1 (c — a)L (c — 0) 

Lemma 4.1. Lei ft e C,ft G C suc/i that Reft > —1, o,kR and e±, 
e 2 >0. 

1. Suppose ei > £2- £ei 

7i(a,&;ft,ft;£i,£ 2 ;s,t) = / d«(-i(u - a) + 2ei) /?1 (-i(u - b) + 2e 2 ) ■ 



Then 

I 1 (a,b;(3 1 ,(3 2 ;e 1 ,e 2 ;s,t) 

(4.6) = Fx (a, 6; ft, ft; ei, £ 2 ; t) - F x (a, b; ft , ft; £1, £ 2 ; s) 
= $i(a,6; ft, ft; £i,£ 2 ;t) - $i(a,6; ft,ft;e 1 ,e 2 ;s), 

with 

Fi(a,b; ft ) ft;ei,e 2 ;t) 

(4.7) = ^-f -fr>> ft+1 (2(£1 - e2 )- i(t - a) )ft 

P2 + 1 

2e 2 -i(t-b) 
2(ei - e 2 ) - z(6 - a) 



x 2 Fx -ft, ft + 1; ft + 2; 



and 



$i(a,6;ft,ft;£i ) e 2 ;i) 

.(2e 2 -i(t-b))^ + fo +l 



(4.8) 



ft +ft + l 

x 2^1 (-ft, -(ft + ft + 1); -ft - ft; 
_ / 2(ei-e 2 ) 



V 2e 2 - i(t - 6) 



2. (Wo hypothesis is needed here.) Let 

ft 

I 2 (a,b;/3i,p 2 ;ei,e 2 ;s,t) = / dii(z(u - a) + 2e 1 ) /3l (-i(u - b) + 2e 2 ) 
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(notice the sign change with respect to Then 
I 2 (a,b; (3i,(3 2 ;e 1 ,e 2 ;s,t) 

(4.9) 

= F 2 (a,b\(3i,f3 2 ;ei,£ 2 ;t) ~ F 2 (a,b;(3 1 ,f3 2 ;ei,s 2 ;s), 

with 

F 2 (a,b;f3i,(3 2 ;ei,e 2 ;t) 

(4 ,o) =i (2 a -i( t - t ))^ (2(£i+£2)+i(t _ a))fl 

Pi + 1 

2e 2 - t(f - 6) 
,2(ei+e 2 ) + i(6-a) 

//, furthermore, a = b = and s,t>0, then 
l2(0,0;P 1 ,p 2 ;si,S2;s,t) 

(4.11) 

= § 2 (!3i,f3 2 -Ei,e 2 -t) - $ 2 {/3i,f3 2 ;si,s 2 ;s), 
^ 2 (Pi,P 2 ;si,e 2 ;t) 

■i ^ 2 -itf^ 

x 2 F, (-p 1 ,-(p 1 + (3 2 + 1); -/3i - /3 2 ; ( ^t^ ) ) ■ 



x 2 Fi -/?i,/3 2 + l;/?2 + 2; 



(4.12) 



Remark. Note the phase factor e l7r/31 in (4.12) is absent from (4.8). It 
is responsible for the divergence of Levy's area for a < 1 /4, as will appear 
clearly during the course of the proof of Theorem 4.4 below. 



Proof of Lemma 4.1. 
1. Set v = u — b. Then, using a deformation of contour, 

h=F 1 {t)-F 1 (s) 

with 

rt-b 

(4.13) Fi{t)= / dv{-iv + 2e 2 f 2 (-i(v + b-a) + 2e 1 f 1 . 

Note that the contour of integration has been implicitly chosen to be a 
line, though any contour going from —2ie 2 to t — b while staying above 
the horizontal line Imz = — 2e 2 would do; this condition is necessary to 
avoid problems of multi-valuedness. Note also that Fi(t) may not be 
defined if (3 2 < — 1 or if e\ < e 2 [the integral may diverge, and the function 
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v — > (i(v + b — a) + 2e\)^ 1 may not be well defined due to problems ofmulti- 
valuedness] . With all this made clear, one may now set x := t Z~b+ 2 i £2 ano - 
obtain 

F 1 (t) = (t-b + 2ie 2 ) 



(4.14) 



x / 1 dx(2e 2 - i(t - b)f 2 x^{2{e l - e 2 ) - i(b - a)) A 



, ( 2g2 -j(t-b))fe+i 

1 p^T[ ( 2 ( £l ~ £2 ) ~ l ( b ~ °)) 

2e 2 - i(t - b) 



(4.15) x 2 Fx -/?i,/3 2 + l;/3 2 + 2; 



2(ei - £2) - »(& - a) 



by formula (4.2). 

Now, by the connection formula (4.3), one gets 

fiM- ^ + ^W .-g) -«»-.))» 

/% + 1 / 2e 2 - i(i - 6) 



(4.16) 



(3 1 +(3 2 + l V2(ei - e 2 ) -i(b-a 
x 2F1 -(1 + ft + ft); -/3i - ft; 

/2(ei -e 2 ) — i(b — a 



+ 



V 2e 2 - i(t - 6) 
r(ft + 2)r(-(/?i+/3 2 + l)) 

r(-/?i) 

2e 2 -i(t-&) \ 



,2(ei - e 2 ) - i(6 - a) 
x 2 F 1 (f3 2 + l,0;2+P 1 +p 2 ; 



(2{e 1 -e 2 )-i{ 



V 2e 2 - i(t - b) 
The second term of the sum simplifies to a constant 
Ci(a,6;/3i,/3 2 ;ei,e 2 ) 

= r(ft + i)r(-(A + A + i)) (2(ei _ _ _ a))A+&+1 

By setting $i(a,6;/?i,/3 2 ;£i,e 2 ;s,t) = Fi(a,6;/3i,/? 2 ;ei,£ 2 ;s,t) - Ci(a,6; 
Pi, fi 2 ; £i) £2)1 one gets the second formula (4.6). 
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2. One may use the same deformation of contour, this time without any 
restriction on the sign of £1 — 62, since Ke(i(v + b — a) + 2e±) > in any 
case. Hence, I 2 (s,t) = F 2 (t) — F 2 (s), with 



(4.17) 



F 2 (t) = (t -b + 2ie 2 ) f 1 dx(2e 2 - i(t - b)) 
Jo 

xx^(2(e 1 +e 2 ) + i(b-a)f 1 
2e 2 - i(t - b) 



x 1 



2(ei +£2) +i(b-a) 



(4.18) 



and, by (4.2), 

F 2 (a,b;0 1 ,0 2 ;e 1 ,e 2 ;t) 

.(2e 2 -t(t-6))^+ 1 



/?2 + l 

x 2 F 1 (-0 1 ,0 2 + l;0 2 + 2; 



(2(e 1+ e 2 ) + i(b-a)f' 
2e 2 - i(t - b) 



2(ei +e 2 ) + i{b-a) 



If one wishes to apply formula (4.3) now as above [see (4.16)], one 
needs to be very careful this time. Calling 



z(t) 



2e 2 - i(t - b) 



k 2(e 1 +e 2 ) + i(6-o) 

the argument of the hyper geometric function, it is quite possible that 
z(t) € M + , in which case the connection formula is not defined. However, 
if a = b = and s,t > 0, then z(t),z(s) G n~ and 

(-zit)) 131 = (z^f 1 • e i7r/ \ (-^(s))^ 1 = (zis)) 131 ■ e^\ 

A computation along the same lines as in the case of 1\ yields 

F 2 (0,0;/3i,/3 2 ;ei,e 2 ;t) 



(4.19) 



C 2 (/3i,/3 2 ;ei,e 2 ) + ie 



^ (2 £2 - 



/3i+/3 2 + l 



x 2^1 -(A + 02 + 1); -A - 021 ( 
where C 2 (0i, 2 ; s\, e 2 ) does not depend on t. 



2(ei+e 2 ; 



2e 2 



□ 



Definition 4.2 (Second moment of Levy's area). Let 



(4.20) 



V(e u e 2 ) 



E 
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If £l = e 2 =: e, then V(ei,e 3 )t = E[(Area 0ii (r(e))) 2 ]. 

The more general quantity considered in Definition 4.2 will be needed 
below for the proof of the convergence of T(e) when e — > in the sense of 
the rough paths. 

Definition 4.3. Let ei,e 2 > 0. We say that ei,£2 — > jointly if £i,£2 — > 
and £i/e2 — ► 1. 

Theorem 4.4. Suppose £i,£2 — >0 jointly in the sense of Definition 
4-3, and a> \ . Then 



(4.21) 

where 

(4.22) C Q 



a(2a-l) 



V(£i,£ 2 )t — >C Q i 



2r(2a-l)r(2a + l) 
r(4a + l) 



la 



+ 



(2a-l)(4a-l) 



Proof. One has 



V(ei,e 2 )t = 2ReE 



(4.23) 



jj'{e l )t^d Xl f\'{e 2 )i^dx 2 

rt rx\ 

xfrr'(£i)-'«%i rr'( £2 ) y + / 2) )^ 



where 



+ 2ReE 



o jo 

" (1 ~ 2a) ") 2 • 2 Re(Vi (£i , e 2 )t + V 2 (ei , e 2 ) t ) 
2 cos 7ra / 



rt rt 2a 

V\ = I dx\ / dyi(-i(x\ - y\) + 2eiY a ~ 
Jo Jo 

r^i ryi 
x / dx 2 / dy 2 (-i(x 2 - y 2 ) + 2e 2 ) 
Jo Jo 



2a-2 



(4.24) 



1 



2a(2a-l)y Jo 



dxi / - y x ) + 2e\) 



2a— 2 



x ((-t(x! - yi ) + 2e 2 ) 2a - (-ixi + 2e 2 ) 



\2a 



(% + 2e 2 ) 2a + (2e 2 



\2a^ 



and 
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N 2a-2 



2a-2 



V 2 = / dxi [ dyx(-i{x x -y{) + 2ei) 
Jo Jo 

x / dx 2 / dy 2 (i(x 2 -2/2) + 2e 2 ) 
Jo Jo 

(4.25) = * _ jf dxi jf - + 2ei) 

x ((i( Xl - yi ) + 2e 2 ) 2a -(ix l + 2e 2 ) 2a 
-(-iy 1 + 2e 2 ) 2a + (2e 2 ) 2a ). 

Computation ofVi- From the last expression (4.24), one gets 
where 

Vi,i = / ctei / (iyi(-i(xi -yi) +2ei) a " {-i{x\ - y\) + 2e 2 ) 
jo Jo 

(4.27) = y* dx(-ia; + 2ei) 2a - 2 (-ix + 2e 2 ) 2Q (t - \x\) 
= 2 Re(tli (0, 0; 2a - 2, 2a; ei , e 2 ; 0, t) - Vi,i, 2 ), 



where 



2a 



Vii2= / dxx(-^ + 2ei) 2Q - 2 (-ix + 2e 2 ) 2a 
Jo 

-|(-it + 2e 1 ) 2Q - 1 t(-it + 2e 2 ) 2a 

- 7i(0,0;2a - l,2a;ei,e 2 ;0,t) 

+ 2ia [ dx((x + £(ei + e 2 )) - + £2)) 
Jo 

x ((-^ + 2ei)(-^ + 2e 2 )) 2a ~ 1 

T |(-it + 2e 1 ) 2Q - 1 t(-^ + 2e 2 ) 2a 

- ii(0,0;2a - 1, 2a; e u e 2 ; 0, t) 

- '-[((-it + 2e 1 )(-it + 2 £2 )) 2a - (4e l£2 ) 2a ] 

+ 2a(ei +e 2 )/i(0,0;2a - 1, 2a - 1; e u e 2 ; 0, t)\. 



2a 
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Suppose s±,£2 — ► jointly (in the sense of Definition 4.3). Then [by using 
the expression (4.6)] 

Re h (0, 0; 2a - 2, 2a; s\ , e 2 ; 0, t) 

l — Rei |(2e 2 - if) 40-1 ^ ^2 - 2a, 1 - 4a; 2 - 4a; 



4a 

(4.28) 



2{ei-e 2 ] 
2e 2 - it 



■ Re i(-ii) 40 - 1 = 3— cos 2^at 4a " 1 . 



4a - 1 4a - 1 

[Note that Im/i(0,0;2a — 2, 2a; e±, e 2 ; 0, i) diverges in the same limit!] 
Next, by (4.6), 

Ji(0,0;2a-l,2a-l;e 1 ,e 2 ;0,t) 

(4.29) x 2 Fi (l - 2a, 1 - 4a; 2 - 4a; ~ ~ ^ ) 

- (2e 2 ) 4a " 1 2i ? i - 2a, 1 - 4a; 2 - 4a; 

hence, (ei + e 2 )/i(0, 0; 2a — 1, 2a — 1; £1, £ 2 ; 0, t) — > 0. Note that this expres- 
sion has no limit if e± and e 2 go independently to 0, and not jointly (in the 
sense of Definition 4.3). 
Finally, 

/i(0,0;2a-l,2a;e 1) e 2 ;0 l t) 

_i(2e 2 -it) Aa ^ (, _. (2{e l -e 2 ) 

(4.30) 



^ l-2a,-4a;l-4a;- , 
4a V V 2e 2 - it 

i(2e 2 ) 4a „ A „ „ , , /2(ei-e 2 ) 



- 2 j 



Fi 1 - 2a, -4a; 1 - 4a; 



4a V V 2e 2 

4a 4a 



Hence, 



(4.31) Vn2^-^-e- 2ma t Aa 

4a 
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and 

(A -to\ V cos27ra 4q 

(432) Vl ' 1 ^-2a(4a-l) t ' 

Let us now turn to the computation of Vi )2 : 

Vi. 2 = f dxx f t dy l (-i(x 1 -y l ) + 2e l ) 2a - 2 (-ix 1 +2e 2 ) 2a 
Jo Jo 

= --^- f dx l {-ix l + 2e 2 ) 2a 
la — 1 Jo 

(4.33) x [(-i(ari - 1) + 2e l ) 2a ~ 1 - (-ia* + 2e 1 ) 2a ^ 1 ] 



■(J 1 (t,0;2a-l,2a;ei,e 2 ;0,t) 



2a- 1 

- ii(0,0;2a - l,2a;ei,e 2 ;0,t)). 

Using (4.6) this time, 

0;2a — l,2a;£i,£ 2 ;0,t) 

* T {(2 £2 - it) 2 - +1 (2(^ - £2 ) + it) 2 "" 1 

2e 2 - it 



2a 

x 2 Fi^l-2a,l + 2a;2 + 2a; 

(4.34) - (2e 2 ) 2a+1 (2(ei - e 2 ) + it) 2 "- 1 

x 2 F 1 ( l-2a,l + 2a;2 + 2a; 



2(ei-e 2 ) + zi 
2e 2 



2(ei -e 2 ) + it 
(-^) 2a+1 (it) 2Q - 1 2 Fi(l - 2a, 1 + 2a; 2 + 2a; 1) 



2a + 1 
i r(2a + l) 2 t4a 
2a T(4a + 1) 

by (4.5). 

The function /i(0, 0; 2a — 1, 2a; e±, e 2 ; 0, t) has already been studied; 
see (4.30). Hence, 

V _ 1 f r(2a)r(2a + l) e~ 2 *™\ 4a 
(435) Vl ^"2a^ll r(4a + l) + "taT ' 



Finally, Vi,3 = Vi, 2 and, clearly, 

vm = ( 4?^iE[r( ei )^( £l )n) (2, 2 ) 2q £1 ^° 4?^?(2^) 2a 



a(l-2a) v " v ;t 7 y a(l-2a) 2 

0. 
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Hence, 



Vi(ei,e 2 ) 



(4.36) 



1 

2a(2a - 1) 

"2r(2a - l)r(2a + 1) 



+ 



COS 27TQ 



ia 



T(4a + 1) (2a-l)(4a-l). 

One may check by using the asymptotic expansion of the Gamma function 
near zero, 

F(e) ~ £ _»o-+7 + 0(£) 
(where 7 is Euler's constant), that this expression is regular when a — > \ or 



a 



2" 



Computation of V 2 ■ The computations are very similar to the previous 
ones, but the final result is surprisingly different, owing to the presence of 
unfortunate (but unevitable) phase factors coming from point 2 of Lemma 
4.1. The reader may easily check step by step that (keeping the same nota- 
tion) 

1 



V 2 



2a(2a 



^y(V 2 ,i-V2,2-V 2 ,3 + V 2 ,4); 



(4.37) 
where 



V 2 ,i = 2 Re(t/ 2 (0, 0; 2a - 2, 2a; e 1 ,e 2 ;0, t) - V 2 ,i, 2 ), 



V 2 ,i, 2 = / dxx{-ix + 2e 1 ) 2a - 2 (ix + 2e 2 ) 2c 
Jo 

= 2^Ti {(-** + tei? a ~ l W + ^2) 



2a 



- I 2 (0,0;2a-l,2a;e 1 ,e 2 ;0,t) 

-2ia [ dx((x + i(e\ - e 2 )) - i(si - e 2 )) 
Jo 

x ((-ix + 2e 1 ){ix + 2e 2 )) 2a - 1 

-{{-it + 2e 1 ) 2a -H(it + 2e 2 ) 2a 

- I 2 (0,0;2a - l,2a;ei,e 2 ;0,t) 

- i[((-tt + 2e x )(it + 2e 2 )) 2a - (4 £l e 2 ) 2a ] 

- 2a{e 1 - e 2 )I 2 (0, 0; 2a -1,2a- 1; e u e 2 ; 0, i) J . 
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Suppose £i,£2 — ► jointly. Then [by using expression (4.11)] 
I 2 (0,0;2a-2,2a;ei,e 2 ;0,t) 
1 



' 2(£i+£ 2 ) 

2e 2 - ^ 



(4.38) x 2-Fi 2 - 2a, 1 - 4a; 2 - 4a; 



-ie 2ina {2e 2 ) ia - 1 



4a- 1 

'2(ei + e 2 ; 



x lim 2 Fi 2 - 2a, 1 - 4a; 2 - 4a; 



2e 2 — is 



The first term in the right-hand side goes to 

1 



t 4a-l 



4a- 1 

as £i, e 2 — ► 0, but the second term is equivalent as e\, e 2 — > while £i/e 2 — >• 1 
to 

(4.39) - — — ie 2i ™(2e 2 ) 4a - 1 lim 2 F X (2 - 2a, 1 - 4a; 2 - 4a; 2 + is), 

whose real part does not miraculously vanish this time, because of the phase 
factor e 2l7ra . However, if a > 4, then this goes to zero. 

The other terms prove to be regular for any a > 0: namely, 

(4.40) (ei -e 2 )/ 2 (0,0;2a- 1,2a- l;e lj e 2 ;0,t) — ► 0; 
(4.41) 



/ 2 (0,0;2a- l,2a;ei,e 2 ;0,t) -c ie i7r(2a_1) — 



4a 



t 4a 

4a 



by (4.12). 
Hence, 

(4.42) v 2 ,i, 2 ^--3-t 4a 



4a 



and, supposing a> \ 



i' 



Let us now turn to the computation of V 2 . 2 : 



) 



(4.44) 



V 2)2 = - ^ T (h (t , 0; 2a - 1 , 2a; e a , e 2 ; 0, t) 

2a — 1 



I 2 (0,0;2a- l,2a;ei,e 2 ;0,t)); 
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using (4.9) this time, 

I 2 (t,0;2a - l,2a;£i,£ 2 ;0,t) 

= 2^Tl{ (2e2 - u ? a+1 ^ + e2 ) - ^ 2a_1 

x 2 Fi 1 - 2a, 1 + 2a; 2 + 2a; 

V V2(ei + £ 2 )-^ 

(4.45) - { 2e 2 ) 2a+1 (2(e 1 +e 2 )-it) 2a - 1 



x 2 F 1 l-2a,l + 2a;2 + 2a; 



2e 2 



2(ei +e 2 ) - it 
(-it) 4a 2 F!(l - 2a, 1 + 2a; 2 + 2a; 1) 



2a + 1 

i r(2o + i)2 _ 2i7ra 



2a T(4a + r 



by (4.5). 
Hence, 

( 4.46) v 22 -^-^-f r ( 2 °> r(2a : l) e 2 - + -V 4 ° 

V ; ' 2a- IV T(4a + 1) 4a/ 

and, z/ a > 4, 

V 2 _ * 

2a(2a - 1) 

(4.47) 

'2cos2 7 raI^H^±l) + 



. la 



T(4a + 1) (2a-l)(4a-l) 
Hence, if a > j, then (after some simplifications) one gets (4.22). □ 

Remarks. Let us compute the constant C a in some particular cases. 
When a — > 1 , C a — ► 4 • 

When a — > i, C a — > ^ [note that the apparent singularity in formula (4.22) 
vanishes] . 

When a ^ \, C a ~ +00. 

Let us conclude this paragraph by mentioning that the piecewise linear 
approximation studied in [5] leads to the same value for C a (in other words, 
for the Levy area of the limiting process) as proved in a recent work by Bau- 
doin and Coutin (see [3]), although a few integration by parts are required 
to see that. Hence, a natural question arises: do both approximations define 
in the limit the same stochastic integration theory? We have no answer for 
the moment, although we believe this should be true. 
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Also, let At ■= JqBs dBs — B(l) s dBg ds be the usual antisymmetric 
Levy area for FBM. Then an easy integration by parts yields E[.4f ] = (4C a — 
l)t 4a . 

4.2. Convergence of Levy's area for a > ^. We shall prove in this sub- 
section the (^-convergence of the q- variation distance d\, d 2 of the Levy area 
for a > \ and q > — . 

We shall be using the following technical lemmas (see [5] and [13]): 

Proposition 4.5. LetO<s<t. 

1. Let w, v be two V -valued paths with q-bounded variation. Let = s + 
(t — s)-^, k = 0, . . . ,2 n , be the nth dyadic subdivision of [s,t]. Then for 
every k > q/2 — 1, there exists a constant C depending only on q and n 
such that 



di(w,v) S)t < C 



En ^ I ^V+?i 4-n V+rc 4-ri I 
^ 1 l k-V Z k l k-l> Z k ' 



n>l k=l 



2. Let w,v be two V -valued paths with q-bounded variation and bounded 
d 2 -norm. Then for every k> q/2—1, there exists a constant C depending 
only on q and k such that 



2" 



n>l k=l 

( 4 - 4 8) +qE» s E Ki^ - 4-^1* 

\n>l fc=l / 



\n>l fe=l / 



1/2 



Lemma 4.6. Lei 
(4.49) <5(K;e;ai,a 2 ) 



■£(|log 2 e|) 



n K 2 n {e ai 2- na2 ) q l 2 



n=0 



2/q 



(E = entire part) for k, s, a± > 0, with a 2 '■= 2a — ot\ . Then 
(4.50) ^25(K;j- l3 ;a 1 ,a 2 ) < oo, 

i/ 2/?(a — -) > 1 and j3ot\ > 1. 
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Proof. One has 

8{n\£\ 0.1,0*2) 



-E(|log 2 e|) 
n=0 



2/q 



Let 7 = (1 — 02(7/2) log 2; the precise estimates depend on the sign of 7. 
Suppose first that 7 < 0: then 

S(k; s; ai, a 2 ) < C(l + | log 2 e |) 2 ( K+1 )/« . e «i , 

hence, X^>i <5( K ; J j a i> a 2) converges if /3ai > 1. 
On the other hand, if 7 > 0, then 

hence, there exists C > (depending on g and k) such that 
6(k; e;ai,a 2 )<C[(l + | log 2 £ \)- £ ^+^/^ £ ^ 
= C(l + \log 2 £\) 2K / q e 2(a - 1/c!) 
and 



(4.51) 



X;*(«;r^;ai,a 2 ) < C^(l + /31ogj-) W 



q~2p(a-l/q) 
J 1 



J>1 J">1 

which converges if 2/3(a — l/q)>l. □ 

Lemma 4.7. Lei fc = 0, ... ,2 n , 6e i/ie ni/i dyadic partition of [0, 1]. 
5ei, /or e > 77 > and k > 0, 



(4.52) 



6'(K;d;e,T}) 



■n>£(|log 2 T/|) fe=l 



where T(e) d G is t/ie tensor d-iterated integral; see (0.4)- 
Then with probability one, 

E5 / (^;^;r /3 ,(j + i)- /3 )<cx), 

i>i 



ifd/3(a-l/q) > 1. 
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Proof. Let, for any choice of indices i\, . . . 



E 



T\e)^dx d 



(4.53) 



: r'^)^)^ [ X1 F'( V )^dx 2 



dx d r'{ v )^d Xd 



Then one may estimate Id very crudely by 



Id(e, V ;t)<C 



(4.54) 
So 



2~ n r2~ n 

dx / dy\ — i(x — y) + 2r/| 
o Jo 



2a— 2 



< C(2- Zn 7] 



2n2a-2\d 



(4.55) 



E[8'(K\e,T))]<C'( n K 2 n (2- 2 V a ~ 2 ) 9/2 ) 

\n>_E(|log 2 ?7|) / 

/ \ d /l 

= C'[ n K 2 n ^\ 1) 

\n>.E(|log 2 Jj|) / 



roo \ d/q 

'1^2^! / 
nd/q d(a— l/q) 



<c"(i+\io g2V \r d ^n 



hence the result. □ 



Corollary 4.8. Suppose there exists a constant C > such that for 
every s,(£f , e > ij > with \s — t\ > Ce, 



(4.56) E[\\T(e)l t - T^f] < £ \t - s 



2a 



(.2 



for some constants Ci > and couples of exponents (a^a^) such that ai 1 + 



oti 2 = 2a and > 0. 



Then the series Y^,j>i(^U )s,t ~ ^(U + l) _/3 )s,t)i.7 — 1) converges a.s. 
in the d2-norm if 2(3(a — l/q) > 1 and f3cti t i > 1, i = 1, . . . , I . Hence, the 
e-regularized Levy area T(e)1 t converges a.s. in the d2-norm when e^O. 
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PROOF. We shall actually prove that E[X)j>i th(T(j-P) 0i t,r((j + l)~ p )o,t)] 
converges. Namely, by Proposition 4.5, 

df(T(j-^,T(U + ir P ))o,t 

/E(f3log 2 j) \ 

(4-57) < E + E 

\ n=0 n>E(f31og 2 j)/ 
2 11 

k=i 

It is a well-known consequence of the hypercontractivity property of the 
Ornstein-Uhlenbeck process (see [16]) that LP-norms (1 < p < oo) are all 
equivalent on any Wiener chaos Ti n (in this case n = 2). Hence, the first 
term in the right-hand side of the last equation is bounded by 

I [/31og 2 j] 
Cj2cf 4 n K 2 n (2- n ) qa ^ ,2 (j- (3 ) qa ^ /2 

i=l n=0 

(4.58) i 

<C'J2(5(K;j-P; ail ,a l2 )y/ 2 

(see Lemma 4.6). 

On the other hand, the second term is bounded by (S'(k; 2;j~@, (j + 
iyP)) q / 2 (see Lemma 4.7). 

Hence (by Holder's inequality, provided q/2 > 1, which is the case when 
a < |), 



E 



■i>i 



< CE f E S(K;j-^,a h ,a i2 )\ + £ «5'(k; 2;^, (j + 1) _/3 ) < 
i=l \ 7 >i / ;>i 



cc 



(by Lemmas 4.6 and 4.7) if 2/3 (a — 1/g) > 1 and /3a^i > 1, i = 1, . . . , I . 

The convergence of the regularized Levy area follows easily along the same 
lines by completeness of the space of rough paths £l q (V). □ 

We shall need below the following definition of the second variation of a 
function with respect to one variable: 



Definition 4.9. Let f :R n — >R be any function of n variables. Then 
the second variation of / with respect to the kth variable is by definition 
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the following function of (n + 1) variables: 

A 2 (/(xi, . . . , x k _i, . . . , x k+1 , . . . , x n ))(e, 77) 

(4.59) 

- 2/1 xi,...,x fc _i,— — ,Xfc+i,...,x n 

+ • • .,a;j fc _i,77 ) x fe+ i, . ..,!„). 

Theorem 4.10. Lei a > j and q > — . T/ien i/te second iterated inte- 
grals T(e)t — T(e) s , Area S] t(r( £ )) converge with respect to the distance d\, 
respectively, d 2 - 

Proof. The quantity one wants to estimate is the following: 
A(e, 77; t) = E[(Areao )t (r(e)) - Area 0it (r(77))) 2 ] 



(4.60) =E 



T l {rj)$dx 1 rT'(r l )^dx 2 



for e > 77 > 0, say. 

By Corollary 4.8, it suffices to prove that, for i>£, A(e,rj;t) is a finite 
sum of terms bounded by a constant times t 2oi2 £ 2ai with ct\ > and ai + 
«2 = 2a. We shall nickname such a term an (01,02) -term for short. Note 
that an (ai, a2)-term is also necessarily an (a'i, a^-term for every a^a^ 
such that a[ + a' 2 = 2a and < a[ < cx\ (our estimates are not necessarily 
optimal). 

Now 

Area , t (r(e)) - Area , t (r(7 7 )) 

V^Wdan f -r'( v )W)dx 2 



XI """J- V" V-JX'z VUX2 

ft rx± 

+ / (r'( £ )W - r'i^hdx, / r'rf) dx 2 



(4.6i) = / r'( £ )W((r(4 2 ) -r(e)P) - (rfa)g> -rrf))^ 

+ ((r( £ )« - r( £ )«) - (rrf) - rrfwrrf) - r(4 2) ; 
V^jgcfaiccr^w -r( £ )«) - (r^w -rrf)) 

£1 ( £ , 77; t) + £ 2 (e, r/; t) + £? 3 (e, 77; t) 
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is the sum of three terms. Note that E^{e,r];t) is of the same form as 
Ei(e,r];t), except for the permutation I^ 1 ) <-> and e <-> Let us es- 
timate E[£'i(e, 77; t) 2 ] and £[£"2(6, 77; t) 2 ] separately. 

Estimation of E2. The first and second component of T decouple and 
one is left with 

E[E 2 (e,r ] ;t) 2 }=E[((T(e)P - T( V )P) - (r( £ )« - r(r?)«)) 2 ] 
(4.62) xE[(T( V )P-r(rj)^) 2 ] 
< C\e - r]\ 2a t 2a 

by Lemma 1.5, which contributes an (a, a)-term to the estimate of Corollary 
4.8, as explained at the very beginning of this proof. 

Estimation of E x . By rewriting (r(e)i 2) - r(e) 2) ) - (r(r ? )i 2 1 ) - T(r]) { Q 2) ) 

as 

(r( e )i 2) -r(,)i 2 ))-(r(4 2 )-r(4 2 )) 



2Re( / Y ,+ ^(z)dz- / T' + ^(z)dz), 

\J xi+in J iri / 



one gets 

(4.63) E[Ei (e, r/; t) 2 } = Wi (e, 77)* + W 2 (e, 7/)t, 



where 



W 1 (e,r ? ) t = 2Re / ek* f "#1 f d Vl 
Jo Jo Jo 

re-T) 

x / d&(-i(xi - yi) + 2e) 

JO 



2a-2 



X ((-i(xi - yi) + fl + 6 + 2??) 



2q-2 



-(-ixi+ 6+6 + 277^ 
- (iyi + 6 + 6 + 2rj) 2a ~ 2 + (a + £ 2 + 2 V ) 2a ~ 2 ) 

= tt~7^ TT' 2Re / / dx 2 (-^(xi - yi) + 2e) Q ~ 

2a(2o! — 1J Jo Jo 

x [(-f(xi - yi) + 2r?) 2a - 2(-i( Xl - Vl ) + e + r?) 2a 

+ H(xi-yi) + 2 £ ) 2a ] 

(4.64) — 7; — 7T- TT- 2Re / dx i / <&2(-t(a:i-l/i) + 2e) 

2a{2a — 1) Jo Jo 



2a-2 
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x [(-ixi + 27]) 2a - 2(-ix! +e + 7]) 2a + (-fii + 2e) 2a ] 
t ft 

. 9Rp / 

2a(2a - 1) 

x [{% Vl + 2 v ) 2a - 2(i Vl + e + V ) 2a + + 2e) 2a ] 

2a-2 



--2Re / dx\ [ dx 2 (—i(x\ — yi) + 2e) 

'a — 1) Jo Jo 



2a-2 



+ /n 1 — — - • 2 Re / dxi [ dx 2 (-i(x 1 -yi) + 2e) 
2a(2a — 1) Jo Jo 

x [(2 v ) 2a -2(e + rj) 2a + (2e) 2a ] 



and 



W 2 (e,ri) t = 2Re / d Xl F Vl / <%i F V d&(-i( Xl - Vl ) + 2e) 2a ~ 2 
Jo Jo Jo Jo 

x ((i(x 1 -y 1 )+6+6+2r ? ) 2a - 2 

- H'yi + + 6 + 2r/) 2Q - 2 + (Ci + £ 2 + 2r / ) 2a ~ 2 ) 

= tt-tt; TT -2Re / / ete2(-i(zi -yi) + 2e) Q " 

2a(2a - 1) Jo Jo 

x - yi ) + 2r ? ) 2a - 2(1(0?! - yi) + £ + r/) 2 ° 

(4.65) + ( i(xi _ yi ) + 2 e ) 2a ] 

~~ „ 77. -r-2Re dxi I dx 2 {-i{xi-yi) + 2e) 2a ~ 2 

2a{2a — 1) Jo Jo 

x [(ixi + 2ri) 2a - 2{ixi +e + T]) 2a + (ix x + 2e) 2a ] 

~ tt^t; TT -2Re / dx i / <£c 2 (-i(xi + 2e) 2a ~ 2 

2a(2a — 1) Jo Jo 

x [(-% + 2?7) 2q - 2(-iyi + e + V ) 2a + (-% + 2e) 2a ] 

1 /"*/"' 2 2 

2a(2a — 1) Jo Jo 

x [(2, ? ) 2a -2( e + r ? ) 2Q + (2 £ ) 2Q ]. 

One has found 

Wi(e, = Vi(e, e) t - 2Vi (e, ^) + Vi(e, = A 2 (V 1 (e, ^ 

(see Definition 4.9 above) and a similar relation for W 2 (£,r])t m terms of 
the second variation of V 2 (see proof of Theorem 4.4 for the definition of 
Vi and V 2 ). Note that the second variation is at most of the same order as 
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the first variation (i.e., the difference between the values at e and rj) — for 
regular functions, it is even of strictly lower order [this is not necessarily 
true otherwise, take, e.g., f(e)=e^, (3 > small!]. Yet all bounds below are 
obtained by estimating the first variation, which is enough for our purposes. 

Let us write asymptotic expansions for e, rj — > 0, e > r), | — ► 1, e, rj <C t of 
each of the terms found in the computation of V\ and V 2 in Theorem 4.4. 

First, 



Re /1 (0,0; 2a - 2,2a;e,r/;0,t) 



(4.66) = — i— Rei{(-it) Aa ~ l (l + 2i(4a - 1)- + O 



_(2-2a)(l-4a)2( £ -,) 

2 - 4a -it V i 2 



hence, 



I A 2 (t Re Ji(0, 0; 2a - 2, 2a; e, •; 0, t))(e, r?)| < Cr/i 4 "" 1 
and similarly for 

|A 2 (iRe/i(0,0;2a- l,2a;e, •;0,t))(e,? ? )|, 

since the arguments of the hypergeometric functions (see 4.31) keep away 
from 1. Next, 

|A 2 ((e + 7/)/i(0,0;2a- 1,2a- 1; e, •; 0, t))(e, r?)| < G^t 4 "" 1 + rf a ); 

|A 2 (t? (-it + 2 £ ) 2a ~ 1 t(-it + 2r ? ) 2a )(e, t?)| < cV^ 1 

and similarly for | A 2 (r/ — > ((— it + 2e)(— it + 2r/)) 2a )(e, r/)| and 

A 2 (r? - (-it + 2 £ ) 2a ~ 1 t(-it + 2 v ) 2a )(e, rj); 

\A 2 (r,^(e V ) 2a )(s, V )\<Cr, 4a . 

In order to estimate I\(t, 0;2a — l,2a;e, 77; 0, , we apply the connection 
formula (4.4) 



2 J 



Fi 1 - 2a, 1 + 2a; 2 + 2a; 



277 -it 



2(e -r)) + it 



T(2 + 2a)r(2a) „ / 2e 
y - 2 Fi 1 - 2a, 1 + 2a; 1 - 2a; 



(4.67) 



+ 



r(l+4a) ' X V '2(e-?7) + it 

2e \ 2a T(2 + 2a)r(-2a) 



2{e-r])+itJ T(l + 2a)r(l -2a) 

2e 



x 2 Fi l + 4a,l;l + 2a;- 

2(e -rj)+it 



A STOCHASTIC CALCULUS BY ANALYTIC EXTENSION 41 

and find out that the term on the first line of (4.34) has a second variation 
bounded by C(t 4a_1 r/ + t 2a r) 2a ), while the second variation of the term on 
the second line is bounded by Ct 2a_1 ry 2a+1 . 
Finally, Vi 4 is an (a, a)-term. 

On the whole, we have proved that Wi(e,rj)t is a sum of (a±, a 2 ) -terms 
for a\ = \ , a, 2a or a + \ ■ 

The same goes for all terms in W^e,^, except of course for the singular 
term in 1/2(0, 0; 2a — 2, 2a; e, 77; 0, t), namely, 

/(£,„) = -_!—. ie 2i ™t(27/) 4a - 1 
4a — 1 

x lim 2 Fi 2 - 2a, 1 - 4a; 2 - 4a; v " 



2rj — is J ' 



see the second line of (4.38). 

The connection formula (4.3) yields 

2(e + r/) 



lim 2-F1 2 - 2a, 1 - 4a; 2 - 4a; 



2r? — is 
r(2 - 4a) 



(4.68) (5+1) e 

\r(l-4a)r(-2a) V r) J 

x 2-P1 [2 - 2a, 1 + 2a; 2 + 2a; 



2o-2 



-2ina 



e + 7] 



T(l + 2a) f e + V Y a - 1 ln(4a _ 1) 
+ T(2-2a)v t? 



Suppose a> \ : then quite simply 



\A 2 f( £ ,r,)\<Cr,^- L t, 
hence, f{e,r]) is a ( 4a 2 ~ 1 » ^)-term. □ 

4.3. Convergence in the d%- distance for \<a< | . An analogue of Propo- 
sition 4.5 holds for the third chaos of FBM, namely: 

Proposition 4.11 (See [5]). Let w,v be two V -valued paths with q- 
bounded variation and bounded d^-norm. Then for every k > q/3 — 1, there 
exists a constant C depending only on q and k such that 



(d 3 (w, v) s , t ) 9/3 < C J2 n K Yl l w |_,^ " V k^ n J q/3 

n>l k=l 
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/ 2" \ 2/3 

+c E»"Eiw|_ lfl -vL fl i ,/2 

\n>l k=l / 
\n>l fc=l / 

(4-69) + C (Y. nK h\™k-^-vk- 1 ^ q ) ' 

\n>l fc=l / 



1/3 



f 2" \ 2/3 

E»"Eiwi, JJ i" 2 + |vi i , i£ i" 2 ) 

^n>l fc=l 



/ 2" \ V d 

\n>l fc=l / 

\n>l fc=l / 



2/3 



Definition 4.12. Let 

Vol 0li (VL,E2,E3) 

(4.70) := fv'{e l )^dx l I^g dx 2 

Jo Jo 

x rv^gw 



The volume functional generated by the e-approximation of a three-dimensional 
fractional Brownian motion is given by Vol(e,e,e), which is more appropri- 
ately written as Vol(r(e)). 

Theorem 4.13. Let a > | and q > -. Then the third iterated integrals 
Vol s> t(r(e)) converge with respect to the distance d%. 

The rest of the article is dedicated to the proof of this theorem. 

Definition 4.14. Let 

, e 2 , e 3 ) t : = E[Vol 0ii (ei, e 2 , £ 3 ) 2 ] 



(4.7i) = e r r ^ • r ^ r 1 ^ ■ / 

~r Jo jo Jo jo jo 



2:2 

dx 3 
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Jo j 
where <r = (eri, a 2 , cr 3 ) £ {±1} 3 - 

Rewrite Vol 0) t(r(e)) - Vol ,t(r(r/)) as 

[Vol ,i(e,e,e) - Vol ,t(e,e,r/)] + [Vol ,t(e,e,r/) -Vol ,t (er, 7/, 77)] 
+ [Vol ,t(e,?7,?7) - Vol ,t(77, 77,77)]. 

Then 

E[(Vol , t (r( e ))-Volo, t (r(r/))) 2 ] 

(4.72) <C(A 2 (W(e,e,-))(s,r 1 )+A 2 (W(e,;r l ))(e,r I ) 

+ A 2 (W(; V , V ))(e,r ] )). 

Corollary 4.8 may be restated in this case: convergence in the <i 3 -norm 
holds true if the following estimate holds: 

(4.73) n\\T(e)l t - T( v )l t \\ 2 } < £ \t - 

i=l 

for some constants Cj > and couples of exponents (ai 1 ,ai 2 ) such that 
«u + Oj 2 = 2a and > 0. 

In the present case (cf. with the definition given in the course of the 
proof of Theorem 4.10), a function f(t,e,rj) will be called an {a\,a 2 )-term 
(qi > 0, a\ + a 2 = 2a) if 

\f(t,e, V )\<Ce 3a H 3a \ 

So one is left with the problem of giving (ai,«2)-type estimates for the 
above A2-terms with e = f] = (j + for some (so that £,77 — > 
and I — > 1). In the following estimates of W(£i,£2,£ 3 )i, one may (and 
shall!) assume that £i,£2,£3 are of the same order, |£ ~ |£ ~ 1, and that 

tS>£i,£2,£3- We shall sometimes write £ = sup(£i,£2,£3) for short. 
The explicit integration 

\ dx 3 / dy 3 (-ia 3 (x 3 - y 3 ) + 2e 3 ) 
Jo Jo 

(4.74) = ((2£ 3 ) 2a - (-ia 3 x 2 + 2£ 3 ) 2a - (ia 3 y 2 + 2£ 3 ) 2a 

+ {~io 3 {x 2 - y 2 ) + 2£ 3 ) 2a )(2a(2a - l))" 1 

yields W(£i, £2, £ 3 )t = Sj=i VVj(£i, £2, £ 3 )t- We shall estimate these four terms 
separately. 



44 J. UNTERBERGER 

Estimation of Wi . Since 

(Of )2ct 

Wi(ex,s 2 ,s 3 )t = -—-^ — T rV(e 1 ,e 2 )t, 
2a{2a — 1) 

this term goes to when ei, e 2 , £3 — > if a > \ and one easily gets {a\,a 2 )- 
type estimates by re-employing those obtained beforehand for Levy's area. 

Estimation of W 2 ■ Two successive integrations by parts yield 

W 2 (ei,£2,£3)t 

= e 2a(Z- r r ^(-^^ - + &i ) 2a_2 



2a-l 



x / dx 2 {-ia 3 x 2 + 2e 3 ) 2a 
Jo 

x [(-«r 2 (aJ2 - yx) + 2e 2 ) 2a - 1 - (-ia 2 x 2 + 2e 2 ) 2a ~ 1 ] 

E o~iQ~2 
r 2a(2a-l) 3 

x / dy 1 {-ia 1 (t-y 1 ) + 2e 1 ) 2a - 1 
Jo 

(4.75) x / dx 2 (-ia 3 x 2 + 2e 3 ) 2a (-ia 2 (x 2 - yi ) + 2e 2 ) 

Jo 

- f dx 1 (-ia 3 x 1 + 2e 3 ) 2a 
Jo 

x [ dyi(-iai(xi-yi) + 2ei) 2a ~ 1 
Jo 

x (-ia 2 (xi - y{) + 2e 2 ) 2a ~ l 

rt 

- / dxi[(-ia\(x\ - t) + 2ei) 2a_1 
Jo 

_(_ i£JlXl + 2ei ) 2a - 1 ] 

x T 1 dx 2 (-ia 3 x 2 + 2e 3 ) 2a (-ia 2 x 2 + 2e 2 ) 2Q ^ 1 
Jo 

which is the sum of three terms, W 2t i, z = 1,2,3. 

(i) (Estimation of W2 1-) 
By Lemma 4.1, 



ft 

/ dx 2 (-ia 3 x 2 + 2e 3 ) 2a (-io- 2 (x 2 - yi ) + 2e 2 ) 
Jo 



2a- 1 
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Sa 3 ,iS a2i ih(0,yi;2a,2a - l;e 3 ,e 2 ;0,i) 



(4.76) + 6 a3 _i£ CT2 -ih(0, yr, 2a, 2a - 1; e 3 , e 2 ; 0, t) 

+ ^ CT3 -15(72,1^(0,1/1; 2a, 2a - l;e 3 ,e 2 ;0,t) 



+ ^3,1^0-2 -ih(0,yi; 2a, 2a - l;e 3 ,e 2 ;0,t). 



If cr 2 = a 3 = 1 and e 3 > e 2 , then ii(0,yi; 2a, 2a - l;e 3 ,e 2 ;0,i) = -Fi(0,yi; 
2a, 2a — l;e 3 ,e 2 ;t) — Fi(0,yi;2a,2a — l;e 3 ,e 2 ;0) (see Lemma 4.1), with 

Fi(0,yi;2a,2a - l;e 3 ,e 2 ;t) 

(4.77) = , ^->l-y^° {2{E3 _ e2) _ iSir 

( 2e 2 - i(t - vi) 

x 2 Fi -2a, 2a; 2a + 1; 



2(e 3 -e 2 ) - iyi 

The real part of the argument z(t) = — ^e^-^-Mix °^ ^ e n yP er g eome t r i c 
function is negative, hence, -Fi(0, y±; 2a, 2a — l;e 3 ,e 2 ;t) is regular in e 3 ,e 2 ,t. 
In the region % ~ 1, one has \z(t)\ <C 1, hence (by Taylor's formula), the 
associated second variations are less than a constant times 

t 

dyi((t - yi) + sf a ~ 1 

(4.78) "° 

x (((* - yi ) + ef a -\ yi + e) 2a + ((t - yi ) + ef a ( yi + e) 2 "" 1 ), 

which has a well defined limit for e — ► as soon as a > g , and one gets a 
(g, 2a — ^)-bound for this term. 

If, on the contrary, ^ < 1, then \z(t)\ is large and one must apply the 
first connection formula (4.3), 

2 Fi(-2a,2a;2a + l;z(t)) 

= r ( 2a + l)( - W> (^-*(\-yi)) 2a 
V 7 \r 2a r l+4a )\2 {e 3 -e 2 )-iyJ 

(4.79) 

/ 1 
x 2 F\ —2a, —4a; 1 — 4a; 



r(-4a) ( 2(e 3 -e 2 )-iyi \ 2a \ 
r(-2a) v 2e 2 -i(t- yi ) J J ' 

The large factor (2(e 3 — e 2 ) — iyi) 2a in the first term of the r.h.s. is compen- 
sated by the overall factor in F±, hence, the same method as before gives a 
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(g, 2a — |)-bound for this term (a > |). The same goes for the second term 
of the r.h.s. 
Now 

Fi(0,yi;2a,2a - l;e 3 ,e 2 ;0) 

(4.80) =i (2(e 3 - e 2 ) - iyi) 

2a 

x 2 F 1 (-2a,2a;2a + l; 



2(e 3 - e 2 ) - % 



can be estimated in the same way in the region y\ = 0{e) [note that Re z(0) is 
not necessarily positive, but | lmz(0)| is bounded from below]. In the region 
yi S> £ where z(0) — 1, the hypergeometric function must be transformed by 
means of the second connection formula (4.4), 

2 Fi (-2a, 2a; 2a + 1; 



2(e 3 -e 2 ) - iyi. 
r(2a + l) / _ 2e 3 



(4.81) 

2e 3 r(-2a - 1 



2(e 3 - £2) ~ iyi 



+ 



2(£ 3 -£ 2 )-i yi T(-2a)T(2a) 

2£ 3 



x 2 Fi 4a + 1,1; 2a + 2; 



2(e 3 -£2) ~ iyi 



which can be estimated by means of Taylor's formula as before, with the 
same result. 

Similar arguments hold if <r 3 = — 1, o 2 = 1; then I 2 (0, yi; 2a, 2a — 1; e 3 , e 2 ; 0, 
t) = F 2 (0,yi;2a,2a - l;e 3 ,e 2 ;t) - F 2 (0, y\\ 2a, 2a - l;e 3 ,e 2 ;0) (see Lem- 
ma 4.1) with F 2 involving the same hypergeometric function but with a 
different argument, namely, wit) = ^^^+1^ ■ The real part of w(t) does 
not have a constant sign this time but keeps away from the forbidden line 
[l;+oo[ (away from the real axis, actually). Details are left to the reader. 

Note that the case e 3 < e 2 can be settled by the same type of arguments, 
for instance, by using an integration by parts to relate to the previous case. 

(ii) (Estimation of W 2]2 .) 

By Lemma 4.1, 

rt 

/ dy 1 (-ia 1 (x 1 - y ± ) + 2£ 1 ) 2a ~ 1 (-ia 2 (x 1 - y x ) + 2£ 2 ) 2a ~ 1 
Jo 

= 5 <T1) _i<5 (r2; _i/i(0,0; 2a - 1,2a - l;ei,e 2 ; -xi,t — x\) 
(4.82) + 5 au i5 a2t ih{<d, 0; 2a - 1, 2a - 1; e^er, -xi,t- x{) 
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+ S aii i5 a2 _iJ 2 (0,0;2a -1,2a- l;e 1 ,e 2 ; -xi,t- xi) 
+ S ai -1^2,1/2(0,0; 2a -1,2a- 1; £1, £ 2 ; ~x\, t - xi). 
In the case o\ = a? = — 1, one writes (see Lemma 4.1) 
li(0, 0; 2a — 1, 2a — 1; £1, £2; —x\,t — xi) 
(4.83) =$i(0, 0;2a- 1,2a- l;£i,£ 2 ;t-xi) 

-$i(0,0;2a- 1,2a- 1; £1, £ 2 ; -xi), 

where 

$i(0, 0; 2a - 1,2a - l;ei,£ 2 ;z) 

(4.84) 

,(2£ 2 -iz) 4a - 1 / 2(£i-£ 2 ) 



i-^-=r 2 Fi 1 - 2a, 1 - 4a; 2 - 4a; 



4a — 1 V 2£2 — iz 

Since the argument of the hypergeometric function keeps small, this case is 
easily settled by using Taylor's formula at order one as before: the associated 
second variations are less than a constant times 



x 4ct-2 n 



(4.85) 



C dx X [x X + £l) 2Q [(xi + £ 2 ) 4a " 2 + ((t - Xi) + £ 2 

Jo 

+ f dX! (Si + El) 2 "- 1 [(XI + E2) 4 "" 1 + ((t-X 1 )+£ 2 
JO 

Hence, one gets once more a (|,2a — |)-bound for this term. 
Suppose now o\ = 1, 02 = — 1: then 

12(0, 0; 2a — 1, 2a — 1; £1, £2; — xi,i — xi) 

(4.86) = F 2 (0, 0; 2a - 1, 2a - 1; e x , e 2 ; t - x x ) 

-F 2 (0,0;2a- 1,2a- 1; £1, £ 2 ; -xi), 



4a- li 



where 



F 2 (0,0;2a-l,2a-l;£i,£ 2 ;z) 



(4-87) =^ (2£ 2 2 a "_ ^ f a (2 ( £ 1+ £2)) 2 - 1 

„ / 2e 2 — 

x 2 Fi 1 -2a, 2a; 2a + 1; 



2(£i+£ 2 )/' 

If z = 0{e), z = x\ or t — x\, then the argument of the hypergeometric 
function is bounded and stays away from [1; +00 [, and 

(4.88) |(2£ 2 - ^) 2a (2( £l + e 2 )) 2a ~ 1 \ < Ce 4a ~\ 

hence, one gets (for a > \ this time) a ( 4c ^~ 1 , 2a % 1 )-bound. Otherwise the 
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first connection formula (4.3) yields 

2e2 — iz 



2 F 1 [ 1 -2a, 2a; 2a + 1; 
r(l + 2a) 



2(ei+e 2 ) 
r(4a-l) / 2(e 1 +e 2 ) xl ~ 2a 



T(2a)r(4a) V 2e 2 - iz 

(4.89) 

2(ei+e 2 ; 



x 2 Fi 1 - 2a, 1 - 4a; 2 - 4a; 



2e2 — iz 

2a 



T(l-4a)f 2(e 1 +e 2 



T(l-2a)V 2e 2 -iz 

The first term of the r.h.s. is similar to that obtained in the case u\ = a 2 = 
— 1, while the second term is a ( , 2a 3 t " 1 )-term if a > |. 

(hi) (Estimation of W2,3.) 

According to the relative signs of 02,(73, 

f Xl dx 2 (-ia 3 x 2 + 2e 3 ) 2a (-ia 2 x 2 + 2e 2 ) 2a ~ l 
Jo 

may be written (up to conjugacy) as 

$i(0,0;2a,2a - l;e 3 ,e 2 ;xi) - $i(0,0;2a,2a - l;e 3 ,e 2 ;0) 

or 

$2(0, 0; 2a, 2a - 1; e 3 , e 2 ; Xi) - lim $2(0, 0; 2a, 2a - 1; e 3 , e 2 ; u). 

?i— >o,«>o 

The argument of the hyper geometric function appearing in $1 remains 
bounded, with a real part <C 1. The second case (involving $2) is slightly 
more complicated; the argument of the hypergeometric function, z(s) = 

2e -fa ' m ust be evaluated for s = x\ and s = 0. If s > Ce 2 , then \z(s)\ 
is fjounded away from 1 and the hypergeometric function is regular and 
bounded. Otherwise one may use the second connnection formula (4.4) as 
in (4.68). Details are left to the reader. 

Estimation 0/W3. It can be obtained from the previous case by simply 
permuting the x-coordinates with the y-coordinates. 

Estimation of W4 . One computes 

W 4 (£l,£2,£3)t 

= dxi dyi(-iai(xi - yi) + 2ei) 2a ~ 2 
Jo Jo 

f Xl f V1 2a~2 2a 

x / dy 2 (-ia 2 (x 2 - y 2 ) + 2s 2 ) (-ia 3 {x 2 - y 2 ) + 2e 3 ) 
Jo Jo 
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r 

- / dxi{-iai(x\ - t) + 2e\Y 
1 JO 

/ dx 2 / dy 2 (-ia 2 (x 2 - y 2 ) + 2e 2 ) 
jo Jo 



2o 

r-.il 



2a— 2 



x (-icr 3 (x 2 - 2/2) + 2e 3 ) 



2a 



2a- 

r-i'i 



+ „ mi -, / dsi / dy 1 (-iai(xi - yi) + 2e 1 ) 2a ~ 
1 Jo Jo 



(4.90) x f 1 dx 2 (-ia 2 (x 2 - yi ) + 2e 2 ) 2a 2 {-ia 3 (x 2 - yi) + 2e 3 ) 

JO 

= - 2a(2 ^_ 1) / f *=iH"l(*i - *) + 2 £l ) 2Q 

x / dy 2 (-ia 2 (x 1 - y 2 ) + 2e 2 ) 2a ~ 2 (-ia 3 (xi - y 2 ) + 2e 3 ) 2c 
Jo 

-M^Ij/o <d2/lHai(t - yi) + 2£l)2a 



x [ dx 2 (-ia 2 (x 2 - yi ) + 2e 2 ) 2a 2 {-ia 3 (x 2 - yi) + 2e 3 ) 2Q 
Jo 

+ ^T7^ tt / ^1/ dyi(—jcri(xi -yi) +2ei) Q 

2a(2a — 1) Jo Jo 

x (-i(7 2 (si - 3/1 ) + 2e 2 ) 2a ~ 2 
x (-ia 3 (xi -yi) + 2e 3 ) 2a . 
The last integral may be rewritten as 

y* dx(-iaix + 2ei) 2Q (-io- 2 x + 2e 2 ) 2Q - 2 (-«T 3 x + 2e 3 ) 2a {t - |x|). 

Since 

([(-ia lX + 2e 1 )(-ia 3 x + 2e 3 )] 2a+1 )' 

= (2a + l)[(-iaix + 2e 1 )(-ia 3 x + 2e 3 )] 2a [-2aia 3 x - 2i(aie 3 + a 3 ex)}, 
one gets 

\2a r 



f dx(-ia lX + 2ei) 2Q (-icj 2 x + 2e 2 ) 2Q - 2 (-icj 3 x + 2e 3 ) 2 

JO 



X 



a\a 3 

ft 



-(crie 3 + (T 3 £i) 



(4.91) x / dx(-ia 1 x + 2e 1 ) 2a (-ia 2 x + 2e 2 ) 2a - 2 (-ia 3 x + 2e 3 ) 2a 

Jo 
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1 

~ 2(2a + l)a 1( 73 

x / dx(-ia 2 x + 2e 2 ) 2a ~ 2 
Jo 

x[((-ia 1 x + 2e 1 )(-ia 3 x + 2e 3 )) 2a+l }'. 
Finally, one is left with the problem of estimating the second variation of 

J 0i0 2 03 : = I dx{-iaix + 2eif 1 {-ia 2 x + 2e 2 ) (3 ' 2 {-ia 3 x + 2e 3 ) (33 
Jo 

with (/?i,/3 2 ,/3 3 ) = (2a,2a-2,2a) or (2o + l,2a-3,2a + l). Taylor's formula 
at order 1 yields very easily a (~, 2a — ^)-bound for a > i. 

REFERENCES 

[1] Abramowitz, M., Stegun, A., Danos, M. and Rafelski, J. (1984). Handbook of 

Mathematical Functions. Harri Deutsch, Frankfurt. MR768931 
[2] Adler, R. J. (1981). The Geometry of Random Fields. Wiley, Chichester. MR611857 
[3] Baudoin, F. and Coutin, L. (2005). Etude en temps petit des solutions d'EDS 

conduites par des mouvements browniens fractionnaires. C. R. Math. Acad. Sci. 

Pans 341 39-42. MR2153390 
[4] Cheridito, P. and Nualart, D. (2005). Stochastic integral of divergence type with 

respect to fractional Brownian motion with Hurst parameter H G (0, j). Ann. 

Inst. H. Pomcare Probab. Statist. 41 1049-1081. MR2172209 
[5] Coutin, L. and Qian, Z. (2002). Stochastic analysis, rough path analysis and 

fractional Brownian motions. Probab. Theory Related Fields 122 108-140. 

MR1883719 

[6] Decreusefond, L. and Ustunel, A. S. (1999). Stochastic analysis of the fractional 

Brownian motion. Potential Anal. 10 177-214. MR1677455 
[7] Dzhaparidze, K. and van Zanten, H. (2004). A series expansion of fractional 

Brownian motion. Probab. Theory Related Fields 130 39-55. MR2092872 
[8] Feyel, D. and de la Pradelle, A. (2001). The FBM Ito's formula through analytic 

continuation. Electron. J. Probab. 6 22. MR1873303 
[9] Gradinaru, M., Nourdin, I., Russo, F. and Vallois, P. (2005). m-order integrals 
and generalized Ito's formula: The case of a fractional Brownian motion with any 
Hurst index. Ann. Inst. H. Pomcare Probab. Statist. 41 781-806. MR2144234 

[10] Kahane, J. -P. (1985). Some Random Series of Functions, 2nd ed. Cambridge Studies 
in Advanced Mathematics 5. Cambridge Univ. Press, Cambridge. MR833073 

[11] Kuhn, T. and Linde, W. (2002). Optimal series representation of fractional Brow- 
nian sheets. Bernoulli 8 669-696. MR1935652 

[12] Lejay, A. (2003). An introduction to rough paths. In Seminaire de Probabilites 
XXXVII. Lecture Notes in Mathematics 1832 1-59. Springer, Berlin. MR2053040 

[13] Ledoux, M., Lyons, T. and Qian, Z. (2002). Levy area of Wiener processes in 
Banach spaces. Ann. Probab. 30 546-578. MR1905851 

[14] Lyons, T. J. (1998). Differential equations driven by rough signals. Rev. Mat. 
Iberoamericana 14 215-310. MR1654527 

[15] Lyons, T. and Qian, Z. (2002). System Control and Rough Paths. Oxford Univ. 
Press, Oxford. MR2036784 



A STOCHASTIC CALCULUS BY ANALYTIC EXTENSION 



51 



[16] Nualart, D. (1995). The Malliavin Calculus and Related Topics. Springer, New 
York. MR1344217 

[17] Pipiras, V. and TAQQU, M. S. (2000). Integration questions related to fractional 
Brownian motion. Probab. Theory Related Fields 118 251-291. MR1790083 

[18] RUSSO, F. and Vallois, P. (1993). Forward, backward and symmetric stochastic 
integration. Probab. Theory Related Fields 97 403-421. MR1245252 

[19] RUSSO, F. and Vallois, P. (2000). Stochastic calculus with respect to continu- 
ous finite quadratic variation processes. Stochastics Stochastics Rep. 70 1-40. 
MR1785063 

[20] van der Vaart, A. W. and Wellner, J. A. (1996). Weak Convergence and Em- 
pirical Processes. Springer, New York. MR1385671 

Institut Elie Cartan 
Laboratoire associe au CNRS UMR 7502 
Universite Henri Poincare Nancy I 
B.P. 239 

F-54506 Vandceuvre les Nancy Cedex 
France 

E-MAIL: jeremie.untcrbcrgcr@iecn.u-nancy.fr 



